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We introduce and study an XY-type model of thermal and quantum phase fluctuations in a two- 
dimensional correlated lattice d-wave superconductor based on the QED3 effective theory of high 
temperature superconductors. General features of and selected results obtained within this model 
were reported earlier in an abbreviated format (Z. Tesanovic, cond-mat/0405235 ). The model is 
geared toward describing not only the long distance but also the intermediate lengthscale physics 
of underdoped cuprates. In particular, we elucidate the dynamical origin and investigate specific 
features of the charge-density-wave of Cooper pairs, which we argue is the state behind the periodic 
charge density modulation discovered in recent scanning-tunneling-microscopy experiments. We 
illustrate how Mott-Hubbard correlations near half-filling suppress superfluid density and favor an 
incompressible state which breaks translational symmetry of the underlying atomic lattice. We show 
how the formation of the Cooper pair charge-density- wave in such a strongly quantum fluctuating 
superconductor can naturally be understood as an Abrikosov-Hofstadter problem in a type-II dual 
superconductor, with the role of the dual magnetic fleld played by the electron density. The resulting 
Abrikosov lattice of dual vortices translates into the periodic modulation of the Bogoliubov-deCennes 
gap function and the electronic density. We numerically study the energetics of various Abrikosov- 
Hofstadter dual vortex arrays and compute their detailed signatures in the single-particle local 
tunneling density of states. A 4x4 checkerboard- type modulation pattern naturally arises as an 
energetically favored ground state at and near the a; = 1/8 doping and produces the local density 
of states in good agreement with experimental observations. 



PACS numbers: 

I. INTRODUCTION 



Several recent experimentsiiSi^ support the proposal 
that the pseudogap state in underdoped cuprates should 
be viewed as a phase-disordered superconductor—. The 
effective theory based on this viewpoint was derived in 
Ref. 13 One starts with the observation^ that in a phase 
fluctuating cuprate superconductor the Cooper pairing 
amplitude is large and robust, resulting in a short co- 
herence length ^ ^ kp^ and small, tight cores for singly 
quantized (anti)vortices. As a result, the phase fluctua- 
tions are greatly enhanced, with hc/2e vortex and an- 
tivortex excitations, their cores containing hardly any 
electrons, quantum tunneling from place to place with 
the greatest of ease, scrambling off-diagonal order in the 
process - incidentally, this is the obvious interpretation 
of the Nernst effect experiments^. Simultaneously, the 
largely inert pairing amplitude takes on a dual respon- 
sibility of suppressing multiply quantized (anti)vortices, 
which require larger cores and cost more kinetic energy, 
while continuously maintaining the pseudogap effect in 
the single electron excitation spectrum. The theory of 
Ref. |5i uses this large d-wave pairing pseudogap A to 
set the stage upon which the low-energy degrees of free- 
dom, identified as electrons organized into Cooper pairs 
and BdC nodal fermions, and fluctuating hc/2e vortex- 
antivortex pairs, mutually interact via two emergent non- 
compact U(l) gauge fields, and a^. v and a couple to 



electronic charge and spin degrees of freedom, respec- 
tively, and mediate interactions which are responsible for 
the three major phases of the theorji^: A c?-wave super- 
conductor, an insulating spin-density wave (SDW, which 
at half-filling turns into a Mott-Hubbard-Neel antiferro- 
magnet) , and an intermediate "algebraic Fermi liquid" , a 
non-Fermi liquid phase characterized by critical, power- 
law correlations of nodal fermions. 

In the context of the above physical picture, the re- 
cent discovery in scanning tunneling microscopy (STM) 
experiments2i2ii£ of the "electron crystal" , manifested 
by a periodic modulation of the local density of 
states (DOS), and the subsequent insightful theoretical 
analysisii of this modulation in terms of the pair density- 
wave, comes not entirely unexpected. Such modulation 
originates from the charge Berry phase term involving 
VQ^i^, the time-like component of v^, and the long- 
distance physics behind it bears some resemblance to 
that of the elementary bosons, like ^He (Ref. IT^ . As the 
quantum phase fluctuations become very strong, they oc- 
casion a suppression of the compressibility of the underly- 
ing electron system, via the phase-particle number uncer- 
tainty relation AipAN > 1, whose effective theory mani- 
festation is precisely the above charge Berry phase. Once 
the off-diagonal order disappears, the system inevitably 
turns incompressible and the diagonal positional order 
sets in, leading to a Mott insulating state. The resulting 
charge-density-wave of Cooper pairs (CPCDW)— causes 
a periodic modulation of the electron density and the size 
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of the pseudogap A and induces a similarly modulated 
local tunneling DOS. In this context, the observed "elec- 
tron crystal" stated should be identified as a CPCDW. 

While the above CPCDW scenario is almost certainly 
qualitatively correct, the ultimate test of the theory is 
whether it can explain and predict some of the specific 
details of the modulation patterns as they are actually 
observed in cuprates. This brings us to the main theme 
of this paper. Typically, when constructing an effective 
low energy theory of a condensed matter system, we are 
solely concerned with the long distance, low energy be- 
havior. In the present case, however, this will not suffice. 
The modulation in question is associated with length and 
energy scales which are intermediate, between the short 
distance scale physics of a single lattice spacing and the 
ultimate long distance behavior. Our aim should thus be 
to construct a description which will be valid not only 
over very long lengthscales but also on the scale of sev- 
eral lattice spacings, which are the periodicities observed 
in experiments2i2ii£. Furthermore, we ideally should be 
aiming for a "bosonized" version of the theory, within 
which a mean-field type approximation for the CPCDW 
state could be gainfully formulated. A natural question 
that needs to be answered first is what should be the 
objects that play the role of these bosons? 

One choice is to designate real-space pairs of elec- 
trons (or holes) as such "elementary" bosons and to en- 
dow them with some "manageable" , i.e. pairwise and 
short-ranged, effective interactions. Such a model indeed 
generically leads to a phase diagram with compressible 
superfluid and incompressible Wigner crystal states just 
as our long-distance argument has suggested. This pic- 
ture of real-space pairs arises within the S0(5)-based 
theorjiMii^, where the low energy sector of the theory 
assumes the form of hard-core plaquette bosons with 
nearest and next-nearest neighbor interactions. When 
one is dealing with extremely strongly bound real space 
s-wave or d-wave pairs this is undoubtedly the natural 
choice. In our view, however, in cuprates one is faced not 
with the real-space pairs but with the momentum space 
Cooper pairs. Apparently, one encounters here an echo 
of the great historical debate on Blatt-Schafroth versus 
BCS pairs - while certain long distance features are the 
same in both limits, many crucial properties are quite 
different^^. To be sure, the Cooper pairs in cuprates 
are not far from the real-space boundary; the coupling is 
strong, the BCS coherence length is short, and the fluctu- 
ations are greatly enhanced. Still, there is a simple litmus 
test that places cuprates squarely on the BCS side: they 
are c?-wave superconductors with nodal fermions. 

This being the case, constructing a theory with Cooper 
pairs as "elementary" bosons turns into a daunting enter- 
prise. Cooper pairs in nodal d-wave superconductors are 
highly non-local objects in the real space, and the effec- 
tive theory in terms of their center-of-mass coordinates 
will reflect this non-locality in an essential way, with com- 
plicated intrinsically multi-body, extended-range interac- 
tions. The basic idea behind the QED3 theory^ is that 



in these circumstances the role of "elementary" bosons 
should be accorded to vortices instead of Cooper pairs. 
Vortices in cuprates, with their small cores, are simple 
real space objects and the effective theory of quantum 
fluctuating vortex-antivortex pairs can be written in the 
form that is local and far simpler to analyze. In this 
dual language the formation of the CDW of Cooper pairs 
translates - via the charge Berry phase discussed above 
- to the familiar Abrikosov-Hofstadter problem in a dual 
superconductoriiSiii. The solution of this problem inti- 
mately reflects the non-local character of Cooper pairs 
and their interactions, and the speciflc CPCDW modu- 
lation patterns that arise in such theory are generally dif- 
ferent from those of a real-space pair density wave. These 
two limits, the Cooper versus the real-space pairs, corre- 
spond to two different regimes of a dual superconductor, 
reminiscent of the strongly type-II versus strongly type-I 
regimes in ordinary superconductors. This difference is 
fundamental and, while both descriptions are legitimate, 
only one has a chance of being relevant for cuprates. 

A well-informed reader will immediately protest that 
the key tenet of the QED3 theory is that we cannot write 
down a useful "bosonized" version of the theory at all - 
the nodal BdG fermions in a fluctuating c?-wave super- 
conductor must be kept as an integral part of the quan- 
tum dynamics in underdoped cuprates. This, while true, 
mostly reflects the central role of nodal fermions in the 
spin channel. In contrast, the formation of CPCDW is 
predominantly a charge sector affair and there, provided 
the theory is reexpressed in terms of vortex-antivortex 
fluctuations - i.e. properly "duahzed" - the effect of 
nodal fermions is less singular and deflnitely treatable. 
This gives one hope that a suitably "bosonized" dual 
version of the charge sector might be devised which will 
provide us with a faithful representation of underdoped 
cuprates. This is the main task we undertake in this 
paper. 

To this end, following a brief review of the QED3 the- 
ory in Section II, we propose in Section III a simple but 
realistic XY-type model of a thermally phase-fluctuating 
d-wave superconductor. Starting from this model we de- 
rive its effective Coulomb gas representation in terms of 
vortex-antivortex pairs. This representation is employed 
in Section IV to construct an effective action for quan- 
tum fluctuations of vortex-antivortex pairs and derive its 
field theory representation in terms of a dual type-II su- 
perconductor, incorporating the effect of nodal fermions. 
In Section V, we discuss some general features of this 
dual type-II superconductor and relate our results to the 
recent STM experiments^'^' This section is written in 
the style that seeks to elucidate basic concepts at the 
expense of overbearing mathematics; we hope the pre- 
sentation can be followed by a general reader. After 
this we plunge in Section VI into a detailed study of 
the dual Abrikosov-Hofstadter-like problem which arises 
in the mean-field approximation applied to a dual super- 
conductor and which regulates various properties of the 
CPCDW state. This particular variant of the Abrikosov- 
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Hofstadter problem arises from the said charge Berry 
phase effect: In the dual representation quantum bosons 
representing fluctuating vortex-antivortex pairs experi- 
ence an overall dual magnetic field, generated by Cooper 
pairs, whose flux per plaquctte of the (dual) Cu02 lattice 
is set by doping: / = p/q = (f — x)/2. Finally, a brief 
summary of our results and conclusions are presented in 
Section VII. 



II. BRIEF REVIEW OF THE QED3 THEORY 

The purpose of this section is mainly pedagogical: Be- 
fore we move on to our main topic, we provide some 
background on the QED3 effective theory of the pairing 
pseudogap in underdoped cuprates. This will serve to 
motivate our interest in constructing an XY-type model 
of fluctuating d-wave superconductors and to make a ca- 
sual reader aware of the connection between the theo- 
retical notions discussed in the rest of this paper and 
the actual physics of real electrons in Cu02 planes. The 
readers well versed in the art of construction of effective 
field theories or those already familiar with the approach 
of Ref. can safely jump directly to Section III. 

The effective theory^ of a strongly fluctuating (1^2 _y2- 
wave superconductor represents the interactions of 
fermions with hc/2e vortex-antivortex excitations in 
terms of two non-compact emergent gauge fields, and 
a^. The Lagrangian C = Cf + Cq is 

(D + iv - ieAY 

^[Dn + ivo - ieAo + ^ - 4- 

Zm 

-iA*'^CT2?7l' + c.c. + /:ob,a] , (1) 

where ^ = ipi), ai's are the Pauli matrices, A is the 
amplitude of the d^2_y2 pairing pseudogap, is the ex- 
ternal electromagnetic gauge field, Z?^ = 9^ + i2a'^^{ai/2) 
is SU(2) covariant derivative and f) = D^ — Dy. Co[v^ a] is 
generated by the Jacobian of the Franz- Tesanovic (FT) 
singular gauge transformation: 

x2-''^5[{wl^)-{nA+nB)]5WI^)-G\B\ , (2) 

where 27m^.s = dx Olpa^b-, w ^ d xv, F{ = ^exf_cuFJ;^^, 
= d^al~d.a^ + 2e^^''a^c4, GXb^ = n{ny^m,, m, = 
(0, 0, (riA — ng)), Q{h) rotates the spin quantization axis 
from z to n, and / dSlfl is the integral over such rotations. 

Its menacing appearance notwithstanding, the physics 
behind ^ is actually quite clear. In which is just 
the effective d-wave pairing Lagrangian, the original elec- 
trons Ca-{x), with the spin quantized along an arbitrary 
direction n, have been turned into topological fermions 
^ = (-01 , ij^i ) through the application of the FT trans- 
formation: (c|,Cj) {exp{i(pA)ip],e'xjp{iipB)ipi), where 



exp(i<y5yi(x) -I- i(pB{x)) equals exp(iip{x)), the center-of- 
mass fluctuating superconducting phase (see Ref. |^ for 
details). The main purpose of this transformation is 
to strip the awkward phase factor exp(iip(x)) from the 
center-of-mass gap function, leaving behind a d-wave 
pseudogap amplitude A and two gauge fields, v and 
a*, which mimic the effects of phase fluctuations. The 
Jacobian given by CQ[v,a^] |2Jl insures that the fluc- 
tuations in continuous fields v and a* faithfully repre- 
sent configurations of discrete (anti)vortices. Note that 
topological fermions do not carry a definite charge and 
are neutral on the average - they do, however, carry 
a definite spin S = ^, reflecting the fact that the spin 
SU(2) symmetry remains intact in a spin-singlet super- 
conductor. As a consequence, the spin density operator 
Sz = C|C| — C|Cj = V'lV'T ~ V'iV'i is an invariant of the 
FT transformation. 

The two gauge fields, v and a', are the main dynami- 
cal agents of the theory. They describe the interactions 
of fermionic BdG quasiparticles with vortex-antivortex 
pair excitations in the fluctuating superconducting phase 
(p{x). A U(l) gauge field is the quantum fluctuating 
superflow and enters the non-conserved charge channel. 
In the presence of the external electromagnetic field A^, 
one has — ^ — eA^ everywhere in >Cy J^l (but not in 
Co), to maintain the local gauge symmetry of Maxwell 
electrodynamics. The Jacobian Cq lj2Jl is the exception 
since it is a purely mathematical object, generated by 
the change of variables from discrete (anti)vortex coordi- 
nates to continuous fields v (and a'). appearing in Cq 
contains only the vortex part of the superflow and is in- 
trinsically gauge invariant. This point is of much impor- 
tance since it is the state of the fluctuating (anti)vortices, 
manifested through Co[v, a'], which determines the state 
of our system as a whole, as we will see momentarily. 

By comparison to v, an SU(2) gauge field aj^ is of 
a more intricate origin. It encodes topological frustra- 
tion experienced by BdG quasiparticles as they encir- 
cle hc/2e vortices; the frustration arises from the ±1 
phase factors picked up by fermions moving through the 
spacetime filled with fluctuating vortex-antivortex pairs. 
These phase factors, unlike those produced by superflow 
and emulated by v, are insensitive to vorticity and de- 
pend only on the spacetime configuration of vortex loops, 
not on their internal orientation: if we picture a fixed 
set of closed loops in 2-f 1-dimensional spacetime we can 
change the orientation of any of the loops any which way 
without affecting the topological phase factors, although 
the ones associated with superflow will change dramat- 
ically. This crucial symmetry of the original problem 
dictates the manner in which a* couples to fermions in 
the above effective theory^. Its topological origin is be- 
trayed by its coupling to a conserved quantum number 
- in spin-singlet superconductors this is a quasiparticle 
spin - and the absence of direct coupling to the non- 
conserved charge channel. By strict rules of the effec- 
tive theory, such minimal coupling to BdG quasiparticles 
is what endows with its ultimate punch at the low- 
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FIG. 1: The phase diagram of cuprates following from the 
QED3 theory. Under the TNomst dome, the Nernst effect 
experiments" indicate strong vortex-antivortex fluctuations; 
this is where the conditions for the applicability of the the- 
ory are met. As the doping is reduced, the superconducting 
ground state is followed by an "algebraic" Fermi liquid, a non- 
Fermi liquid state described by the symmetric phase of QED3 
and characterized by critical, power-law correlations^^ . At yet 
lower doping, the BdG chiral symmetry is broken, resulting in 
an incommensurate antiferromagnet (SDW), which eventually 
morphs into the Neel antiferromagnet at half-filling. Depend- 
ing on material parameters (anisotropy, strength of residual 
interactions, degree of interlayer coupling, etc.) of a particu- 
lar HTS system, the transition between a superconductor and 
a SDW could be a direct one, without the intervening chiral 
symmetric ground state. The way CPCDW, the main subject 
of this paper, fits into this phase diagram is detailed later in 
the text. 



est energies. While in this limit - which couples to 
fermions in the non-minimal, non-gauge invariant way 
- turns massive and ultimately irrelevant, ends up 
generating the interacting infrared (IR) critical point of 
QED3 theory which regulates the low energy fermiology 
of the pseudogap state. 

£ manifestly displays the U(l) charge and SU(2) 
spin global symmetries, and is useful for general consid- 
erations. To perform explicit calculations in the above 
theory one must face up to a bit of algebra that goes 
into computing Co[v,a^]. The main technical obstacle 
is provided by two constraints: i) the sources of 
should be an SU(2) gauge equivalent of half-flux Dirac 
strings, and ii) these sources are permanently confined 
to the sources of , which themselves are also half-flux 
Dirac strings. The confined objects formed by these two 
half-string species are nothing but the physical fluctuat- 
ing hc/2e vortex excitations. These constraints are the 
main source of technical hurdles encountered in trying 
to reduce (^^l to a more manageable form - nevertheless, 
it is imperative they be carefully enforced lest we lose 



the essential physics of the original problem. The con- 
straints can be solved by introducing two dual complex 
fields = l^le*-^ and <i>s{x) = |$,|e*"*= to provide 

coherent functional integral representation of the half- 
strings corresponding to and a^, respectively. In this 
dual language the half-strings are simply the worldlines 
of dual bosons $ and $s propagating through (2 -I- 1)- 
dimensional spacetime. It is important to stress that 
this is just a mathematical tool used to properly cou- 
ple vortex-antivortex pair fluctuations to the electronic 
degrees of freedom - readers less intimate with the tech- 
nology of dualization will find it explained in detail in 
the Appendix. The confinement of the two species of 
half-strings into hc/2e (anti)vortices is accomplished by 
demanding |<I'(2:)| = |<I>s(a;)|. This finally results in £0 
expressed as 

exp{- J^^ dT J d\Co[v^,al]) J dfln J V[^, Ad, k'] 

X C-i[|$|] exp{ J d^x{i2Ad ■ w + i2K' ■ - Cd)) , 

(3) 

where Cd[^, Ad, k*] is a dual Lagrangian given by 

mlm' + \i^,~^2nAd,m' + ^m'' 

+ |<i>P(a^0, - 2^f^^3^;)2 , (4) 

and C[|<I'|] is a normalization factor determined by 

J dnnV[a\(l)s,K']exp{i2K'-F'+\^f{d^(t)s~2T:n'^Klf} . 

(5) 

Here Ad and are the charge and spin dual gauge fields, 
respectively, having been introduced to enforce the two 
i5-functions in Relative to our notation here, in Ref. 
0both gauge fields are rescaled as ^^(k') 2T:Ad{K^)- 
Notice that Ad and couple to v and a' as exp(i2Ad ■ 
w + i2K''-F'-) = exp[i2Ad-{dxv) + iK\exf_„y{df_tal-d^a^^ + 
2e'^^alal)) ®. 

At this point, one should take note of the following 
convenient property of £ H1I3|I . one that will prove handy 
repeatedly as we go along: Ad and provide - via the 
above coupling to v and a' - the only link of commu- 
nication between the fermionic matter Lagrangian £y 
and the vortex-antivortex state Jacobian £9. This way 
of formulating the theory is more than a mere math- 
ematical nicety: It is particularly well-suited for the 
strongly fluctuating superconductor which nonetheless 
does not belong to the extreme strongly bound limit of 
real-space pairs, when they can be viewed as "elemen- 
tary" bosons. Such extreme limit of "preformed" pairs 
or pair "molecules" is frequently invoked, inappropriately 
in our view, to describe the cuprates. Were the cuprates 
in this extreme limit, there would be no gapless nodal 
quasiparticle excitations, the situation manifestly at odds 
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with the available experimental information^^. In addi- 
tion, such demise of nodal fermions would rob the effec- 
tive low-energy theory of any spin degrees of freedom, a 
most undesirable theoretical feature in the context of all 
we know about the cuprates. 

In a superconductor, the vortex-antivortex pairs are 
bound and the dual bosons representing vortex loop 
worldlines are in the "normal" (i.e., non-superfluid state): 
($) = (<i)^) — (see the Appendix). Consequently, both 
dual gauge fields are Ad and are massless; this is as 
it should be in the dual normal state. This implies that 
upon integration over Ad and both and aj^ will 
be massive: Cq M^Vq + M^v^ + {v — + a), where 
Mq oc Ml ^ £^d, the dual correlation length. As ad- 
vertised above, the massive character of v immediately 
makes the system a physical superconductor: If we in- 
troduce a static external source vector potential A*^^ (q) , 
the response to this perturbation is determined by the 
fermionic "mean- field" stiffness originating from Cf and 
the "intrinsic" stiffness set by Co, whichever is smaller. 
In a strongly fluctuating superconductor this in practice 
means Co- The mass M|(x,T) of the Dopplcr gauge 
field V in Cq is effectively the helicity modulus of such 
a superconductor and determines its superfluid density 
Ps, reduced far bellow what would follow from Cf alone. 
In a similar vein, the mass Mq(x,T) of Vq sets the com- 
pressibility of the system Xc, again far smaller than the 
value for the non-interacting system of electrons at the 
same density: 



lim — 
q^o ,5A^^(q)(5A^^(-q) 

~^q™ M-(q)M-(-q) 



Jl]^{q)Ml 



n}^(q) 



Ml 



nf(q)M„2 



nf(q) 



Ml 



Mi 



M' 



(6) 

where Eq is the ground state energy and 11^"' (q; x, T) and 
n^''(q;x,T) are the current-current and density-density 
responses of the fermionic Lagrangian C / , respectively. 

Evidently, the "intrinsic" stiffness of v as it appears 
in Cq encodes at the level of the effective theory the 
strong Mott-Hubbard short-range correlations which are 
the root cause of reduced compressibility and superfluid 
density and enhanced phase fluctuations in underdoped 
cuprates. This is entirely consistent with the spirit of the 
effective theory and should be contrasted with the ap- 
plication of a Gutzwiller-type projector to a mean-field 
BCS state, to which it is clearly superior: In the latter 
approach one is suppressing on-site double occupancy by 
a local projector mimicking correlations that are already 
included at the level of the renormalized "mean-field" 
fermionic Lagrangian Cf in (Q, i.e. they are built-in 
into Ily' and 11^'' ® . In this approach the superconduc- 
tivity is always present at T = as long as the dop- 
ing is finite, albeit with a reduced superfluid density. 
In contrast, within the effective theory (Q, the micro- 
scopic Mott-Hubbard correlations are additionally echoed 
by the appearance of well-defined vortex-antivortex exci- 



tations, whose core size is kept small by a large pairing 
pseudogap A. Such vortex-antivortex excitations further 
suppress superfluid density from its renormalized "mean- 
field" value and ultimately destroy the off-diagonal long 
range order (ODLRO) altogether - even at T = and 
finite doping - all the while remaining sharply defined 
with no perceptible reduction in A. 

As the vortex-antivortex pairs unbind, the supercon- 
ductivity is replaced by the dual superfluid order: ($) ^ 
0, ($5) ^ 0, and and turn massless: £0 
co?(Vwo)^+cC(Vxv)2 + (u ^ a), wheref(a;,T ^ 0) is the 
superconducting correlation length, cq and c are numeri- 
cal constants, and we have used radiation gauge V-v = 0. 
Physically, the massless character of and describes 
the admixing of free quantum vortex-antivortex excita- 
tions into the ground state of the system which started as 
a d-wave superconductor. Now, it is clear from Eq. © 
that the response to A'^^(q) is entirely determined by Cq. 
The massless character of v in Cq implies the vanishing of 
the helicity modulus and superfluid density, despite the 
fact that the contribution from Cf to both remains finite 
and hardly changes through the transition. Similarly, the 
response to ^^^(q) vanishes as well, since vq is now mass- 
less. Thus, the system simultaneously loses superfluidity 
{ps — > 0) and turns incompressible (xc coC|qp 0) 
for doping smaller than some critical Xc- 

Note that the theory^ predicts a universal relation be- 
tween the superfluid density ps and the fluctuation dia- 
magnetic susceptibility Xdm in underdoped cuprates, as 
T — !■ 0. Right before the superconductivity disappears, 
for X > Xc, the discussion surrounding Eq. (|SJ| implies 
Ps ^ Ml ^ Similarly, in the region of strong super- 

conducting fluctuations for x < Xc, Ml in Q is replaced 
by ofq^ as is clear from the previous paragraph. The 
prefactor is nothing but the diamagnetic susceptibility 
Xdia ~ c^- Assuming that the superconducting correla- 
tion length and its dual are proportional to each other, 
i ^ S,d, one finally obtains Xdia ~ Pj^- In the case where 
strong phase fluctuations deviate from "relativistic" be- 
havior, the dynamical critical exponent z ^ \ needs to 
be introduced and the above expressions generalized to 

Ps ^ Ml ^ Xdia ~ and Xdia ~ Ex- 

perimental observation of such a universal relation be- 
tween Ps and Xdia would provide a powerful evidence for 
the dominance of quantum phase fluctuations in the pseu- 
dogap state of underdoped cuprates. 

In the "dual mean-field" approximation (<f>) = ($s) = 
$ one finds: 



Cq 



1 



47r2|$|2 

^-Jd^xc^^(2-^J x)K'exp{J d^x{i2K^ ■ +n[K^])} , 

(7) 



where / d^xTZ = In / dr^fl exp(-47r2 / d^x\^\'^{K 

£g[a'] appears somewhat unwieldy, chiefly through its 
non-local character. The non-locality is the penalty we 
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pay for staying faithful to the underlying physics: While 
the topological origin of a' demands coupling to the con- 
served SU(2) spin three-currents, its sources are perma- 
nently confined to Dirac half-strings of the superflow 
Doppler field v, which by its very definition is a non- 
compact U(l) gauge field (see below). The manifestly 
SU(2) invariant form subjected to these constraints 
is therefore even more elegant than it is useful. 

The situation, however, can be remedied entirely by a 
judicious choice of gauge: = = 0, = a^. In this 
"spin-axial" gauge - which amounts to selecting a fixed 
(but arbitrary) spin quantization axis - the integration 
over Ki can be performed and Cq (|7|) reduces to a simple 
local Maxwellian: 



Co [v,a] 



1 



47r2|$|2 



{d X vf + 



1 



47r2|$|2 



{d X a) 



(8) 



where $ is the dual order parameter of the pseudogap 
state, i.e. the condensate of loops formed by vortex- 
antivortex creation and annihilation processes in (2+1)- 
dimensional Euclidean spacetime. This is the quan- 
tum version of the Kosterlitz-Thouless unbound vortex- 
antivortex pairs and is discussed in detail in the Ap- 
pendix. The effective Lagrangian £ = £f + £o[v,a] of 
the quantum fluctuating d-wave superconductor finally 
takes the form: 

(B + iv-ieA)'^ 

* Do + ivo - ieAo + ^ - ^ * 

2m 

- i A5''^cr2?}1' + c.c. + £0^,0] , (9) 

where — > 9^ + ja^crs, other quantities remain as de- 
fined below Eq. QJ, and £o[v,a] of the pseudogap state 
is given by Eq. More generally, the full spin-axial 

gauge expression for £o[?^7a] in ©7 valid both in the 
pseudogap ((<i>) ^ 0) and superconducting (($) — 0) 
states, is: 

J dr j (fr£o[v^,a^]) J Vl^, (j)^, A^, n] 

X C^^[|$|] exp{ J (fx{i2Ad ■ w + i2K ■ F ~ £d)} , 

(10) 

where dual Lagrangian £d[^, Ad, n] is 

mlm' + \{^,~^2nAd,)<^>\' + ^m' 

+ mHd^^,-27:K,,f , (11) 
and normalization factor C[|<i>|] equals 

J I?[a,(/.s,K]exp{z2K-F-|-|$P(a^0, ~27r«;^)2} . (12) 

Again, Ad and k are the charge and spin dual gauge fields, 
respectively, which, in the spin-axial gauge, couple to v 
and a as exp(i2Ad ■ w + i2K ■ F) = exp(i2Ad ■ {d x v) + 



The above is just the standard form of the QED3 the- 
ory discussed earlier^-. It resurfaces here as a particular 
gauge edition of a more symmetric, but also far more 
cumbersome description - fortunately, in contrast to its 
high symmetry parent, the Lagrangian ^ itself is emi- 
nently treatable^ii. 

A committed reader should note that the ultimate non- 
compact U(l) gauge theory form of £ © arises as a natu- 
ral consequence of the constraints described above. Once 
the sources of v and a are confined into physical hc/2e 
vortices, the non-compact U(l) character is the shared 
fate for both. While we do have the choice of selecting the 
singular gauge in which to represent the Berry gauge field 
a, there is no similar choice for u, which must be a non- 
compact U(l) gauge field. It is then only natural to use 
the FT gauge and represent a as a non-compact U(l) field 
as well, in order to straightforwardly enforce the confine- 
ment of their respective sourceaSi. The non-compactness 
in this context reflects nothing but an elementary prop- 
erty of a phase- fluctuating superconductor: The conser- 
vation of a topological vortex charge d-n = = d-nA,B^- 

We end this section with the following remarks: In ex- 
plicit calculations with (jS)) or its lattice equivalent, it is 
often useful to separate the low energy nodal BdG quasi- 
particle excitations from the rest of the electronic degrees 
of freedom by linearizing £/ near the nodes. The nodal 
fermions 4'u,a, where a = 1, I, 2 and 2 is a node in- 
dex, can then be arranged into N four-component BdG- 
Dirac spinors following the conventions of Ref. |5j, where 
N = 2 for a single CUO2 layer. These massless Dirac-like 
objects carry no overall charge and are at zero chem- 
ical potential - reflecting the fact that pairing in the 
particle-particle channel always pins the d-wave nodes 
to the true chemical potential - and can be thought of 
as the particle-hole excitations of the BCS "vacuum" . 
They, however, can be polarized and their polarization 
will renormalize the fluctuations of the Doppler gauge 
field V, the point which will be emphasized later. Fur- 
thermore, the nodal fermions carry spin S = ^ and in- 
teract strongly with the Berry gauge field a to which they 
are minimally coupled. In contrast, the rest of the elec- 
tronic degrees of freedom, which we label "anti-nodal" 
fermions, '!/'cr,{Q/3) ; where (a/3) = (12), (21), (12), and 
(21), are combined into spin-singlet Cooper pairs and do 
not contribute significantly to the spin channel. On the 
other hand, these anti-nodal fermions have finite density 
and carry the overall electric charge. Their coupling to v 
is the dynamical driving force behind the formation of the 
CPCDW. Meanwhile, the nodal fermions are not affected 
by the CPCDW at the leading order - their low energy 
effective theory is still the symmetric QED3 even though 
the translational symmetry is broken by the CPCDWi^. 

The presence of these massless Dirac-like excitations 
in is at the heart of the QED3 theory of cuprates^. 
QED3 theory is an effective low energy description of a 
fluctuating d^2_y-2-waNe superconductor, gradually los- 
ing phase coherence by progressive admixing of quantum 
vortex-antivortex fluctuations into its ground state. All 
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the while, even as the ODLRO is lost, the amplitude of 
the BdG gap function remains finite and largely undis- 
turbed as the doping is reduced toward half-filling as de- 
picted in Fig. n How is this possible? A nodal d^2_y2- 
wave superconductor - in contrast to its fully gapped 
cousin or a conventional s-wave superconductor - pos- 
sesses two fundamental symmetries in its ground state; 
the familiar one is just the presence of the ODLRO, 
shared by all superconductors. In addition, there is a 
more subtle symmetry of its low energy fermionic spec- 
trum which is exclusively tied to the presence of nodes. 
This symmetry is emergent, in the sense that it is not a 
symmetry of the full microscopic Hamiltonian but only 
of its low energy nodal sector, and is little more than the 
freedom intrinsic to arranging two-component nodal BdG 
spinors into four-component massless Dirac fermions^i - 
by analogy to the field theory we call this symmetry chi- 
ral. The QED3 theory first formulates and then answers 
in precise mathematical terms the following question: 
can a nodal (i3,2_j,2-wave superconductor lose ODLRO 
but nonetheless retain the BdG chiral symmetry of its low 
energy fermionic spectrum^? In conventional BCS theory 
the answer is a straightforward "no" : as the gap goes to 
zero all vestiges of the superconducting state are erased 
and one recovers a Fermi liquid normal state. However, 
in a strongly quantum fluctuating superconductor con- 
sidered here, which loses ODLRO via vortex-antivortex 
unbinding, the answer is a remarkable "yes" . The chirally 
symmetric, IR critical phase of QED3 is the explicit re- 
alization of this new, non-Fermi liquid state of quantum 
matter, with the phase order of a d^^-y^-vave supercon- 
ductor gone but the chiral symmetry of fermionic exci- 
tations left in its wake. While these excitations are in- 
coherent, being strongly scattered by the massless Berry 
gauge field, the BdG chiral symmetry of the low energy 
sector remains intacliS^. Importantly, as interactions get 
too strong ~ the case in point being the cuprates at very 
low dopings - the BdG chiral symmetry will be eventually 
spontaneously broken, leading to antiferromagnetism and 
possibly other fully gapped states^'^. Nonetheless, the 
BdG chiral symmetry breaking is not fundamentally tied 
to the loss of ODLRO: different HTS materials are likely 
to have different T — Q phase diagrams with larger or 
smaller portions of a stable critical "normal" state - de- 
scribed by the symmetric QED3 - sandwiched between a 
superconductor and an antiferromagnet (SDW) (Fig. ^1. 
In all cases, provided our starting assumption of the pre- 
dominantly pairing nature of the pseudogap is correct, 
the symmetric QED3 emerges as the underlying effective 
theory of underdoped cuprates, echoing the role played 
by Landau Fermi liquid theory in conventional metals. 



III. THERMAL ("CLASSICAL") PHASE 
FLUCTUATIONS 

Our first goal is to introduce an XY model type repre- 
sentation of thermal (or "classical") phase fluctuations. 



To this end we first observe that in cuprates, we are deal- 
ing with a d-wave superconductor on a tightly bound 
two-dimensional Cu02 lattice (the black lattice in Fig. 
121. The simplest starting point is the lattice d-wave su- 
perconductor (LdSC) model discussed in detail by Vafek 
et al}^-. The model describes fermions hopping between 
nearest neighbor sites {ij) on a square lattice with a 
renormalized matrix element t* and contains a nearest 
neighbor spin-singlet pairing term with an effective cou- 
pling constant Aofi adjusted to stabilize the (i^2_j,2 state 
with the maximum pairing gap A: 

-ffLdSC = -t* c\„Cja + ^ Ay [cj^c]^ - c|^c]-f] + (C.C.) 

-f(l/Acff)5]|A,,f -^(...) , (13) 

where (• • • ) denotes various residual interaction terms. 
c|g. and Cia are creation and annihilation operators of 
some effective electron states, appropriate for energy 
scales below and around A, which already include renor- 
malizations generated by integration of higher energy 
configurations, particularly those associated with strong 
Mott-Hubbard correlations. This effective LdSC model is 
phenomenological but can be justified within a more mi- 
croscopic approach, an example being the one based on a 
t — J-style effective Hamiltonian. The complex gap func- 
tion is defined on the bonds of the Cu02 (black) lattice: 
Aij — Aij e:xjp(i9ij) . The amplitude Aij is frozen be- 
low the pseudogap energy scale A and equals Aij = ±A 
along horizontal (vertical) bonds. Thus, the d-wave char- 
acter of the pairing has been incorporated directly into 
Aij from the start. 

What remains are the fluctuations of the bond phase 
6ij. It is advantageous to represent these bond phase 
fluctuations in terms of site fluctuations by identifying 
exp(i6'y) e:xp(iipk), where (pk is a site phase variable 
associated with the middle of the bond {ij). Conse- 
quently, we are now dealing with the set of site phases 
(fik located at the vertices of the blue lattice in Fig. |21 
Note that the blue lattice has twice as many sites as the 
original Cu02 lattice or its dual. This is an important 
point and will be discussed shortly. 

We can now integrate over the fermions in the LdSC 
model and generate various couplings among phase fac- 
tors exp{i(pk) residing on different "blue" sites. The re- 
sult is the minimal site XY-type model Hamiltonian rep- 
resenting a fluctuating classical lattice d-wave supercon- 
ductor: 

nn rnnn 

- J2 ^ cos{^, - ipj) , (14) 

hnnn 

where X^nn runs over nearest neighbors on the blue lat- 
tice while X^rnnn ^^'^ Sbnnn '^^'^ ^'^'^r red and black 
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FIG. 2: Cu02 lattice represented by thick black lines (ver- 
tices are copper atoms and oxygens are in the middle of each 
bond). The nodes of the blue lattice (thin solid lines) repre- 
sent sites of our "effective" XY model l|14|l . with the phase 
factor exp(i</5i) located at each blue vertex. The dual of the 
copper lattice is shown in red dashed lines. The nearest neigh- 
bor coupling on the blue lattice is denoted by J (blue link) 
and the two next-nearest neighbor couplings are Ji (red link) 
and J2 (black link) . Generically, Ji 7^ J2 and the translational 
symmetry of the blue lattice is broken by the "checkerboard" 
array of "red" and "black" plaquettes, as is evident from the 
figure. This doubles up the unit cell which then coincides 
with that of the original (black) CUO2 lattice, as it should be. 



next-nearest neighbor bonds depicted in Fig. |21 respec- 
tively. The couplings beyond next-nearest neighbors are 
neglected^i in ifT^ . Such more distant couplings are not 
important for the basic physics which is of interest to 
us here, and their effect can be simulated by a judicious 
choice of Ji and J2 - the reader should be warned that 
the situation changes when T — > as will be discussed in 
the next section. Keeping Ji and J2, however, is essen- 
tial. If only J is kept in H14|l . the translational symmetry 
of the blue lattice would be left intact and we would find 
ourselves in a situation more symmetric than warranted 
by the microscopic physics of a d-wave superconductor. 
In this sense Ji and J2 or, more precisely, their differ- 
ence Ji — J2 ^ 0, are dangerously irrelevant couplings 
and should be included in any approach which aims to 
describe the physics at intermediate lengthscales. 

In general, near half-filling, J, Ji and J2 tend to be 
all positive and mutually different. This ensures that 
the ground state is indeed a d-wave superconductor (the 
reader should recall that d-wave signs had already been 
absorbed into Ay's so the ground state is an XY ferro- 
magnet instead of an antiferromagnet - the antiferromag- 



net would be an extended s-wave state.). The explicit 
values of J, Ji and J2 can be computed in a particular 
model. In this paper we treat them as adjustable param- 
eters. 



A. Vortex-antivortex Coulomb gas representation 

of 

The fact that Ji ^ J2 breaks translational symmetry 
of the blue lattice. This is as it should be since the blue 
lattice has twice as many sites as the original Cu02 lat- 
tice. This is an important point to which we will return 
shortly. Now, however, let us assume for the moment that 
J I = J2. This assumption is used strictly for pedagogical 
purposes. Then, the blue lattice has a full translational 
symmetry and it is straightforward to derive the effec- 
tive Coulomb gas representation for fluctuating vortex- 
antivortex pairs. Here we follow the derivation of Ref. 
I25I The cosine functions in ((TH) are expanded to second 
order and one obtains the effective continuum theory: 

■Hcout^\j j d?r\V'p{v)\^ + {■■■) , (15) 

where J = J + Ji + J2 and ^p{r) is the continuum version 
of the site phase on the blue lattice. (• • • ) denotes 
higher order terms in the expansion of the cosine func- 
tion. As usual, the superfluid velocity part u(r) — V(y3(r) 
is separated into a regular (XY spin-wave) part and a 
singular (vortex) part (V^s)!,: 

t;(r) = Vx(r) + 2tt{z j d^r'n{r')G{r, r') , (16) 

where x is a regular free field and n(r) is the density of 
topological vortex charge: n(r) =J2i^{^~^i)~J2jH'''^~ 
Tj), with and being the positions of vortex and an- 
tivortex defects, respectively. We are limiting ourselves 
to ±1 vortex charges but higher charges are easily in- 
cluded. G{r,r') is the 2D electrostatic Green's function 
which satisfies: 

V2G(r,r') = 6{r-r') . (17) 
Far from systems boundaries, the solution of (|17|l is: 

Gir,r')^^H\r-r'\/a) + C , (18) 

where a is the UV cutoff, of the order of the (blue) lattice 
spacing and C is an integration constant, to be associated 
with the core energy. 

This ultimately gives the vortex part of the Hamilto- 
nian as: 

Hy = 27r2j j(fr(fr'(fr"n{Y')n{Y")VG{r,r')-VG{r,r") , 

(19) 
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After integration by parts this results in the desired 
Coulomb gas representation of Ti.^ : 

-IT J j (f r(fr'n{r)n{r') ln{\r - r'\/ a) + Ec J (frn^{r) , 

(20) 

where Ec ^ C is the core energy. Note that the integrals 
in the first part of (|20|) include only the regions outside 
vortex cores (the size of which is ~ a^). 

How is the above derivation affected if we now restore 
Ji 7^ J2, as is the case in real cuprates? The long dis- 
tance part remains the same since Ji — J2 enters only 
at the 0{k'^) order of the gradient expansion of cosine 
functions in H14|) . i.e. Ji — J2 7^ affects only the terms 
denoted by (• • • ) in Tlcont CSI- So, the strength of the 
Coulomb interaction between (anti)vortices in TCv H2U|) is 
still given by J = J + Ji + J2. The core energy changes, 
however. This change can be traced back to (|17|l and 
()18|l : depending on whether the vortex core is placed in 
a red (i.e., containing the red cross in Fig. [21) or a black 
plaquette of the blue lattice, the constant of integration 
C will generally be different. This is just the reflection of 
the fact that for Ji ^ J2 the original translational sym- 
metry of the blue lattice is broken down to the checker- 
board pattern as is obvious from Fig. |2 So, Eq. H18() 
must be replaced by: 



G(r,r')-;^ln(|r- 



c'\/a) + Cr 



(b) 



(21) 



G(r, r') is just the electrostatic potential at point r pro- 
duced by a vortex charge at r'. Since not all locations for 
r' are equivalent, there are two constants of integration 
Cr 7^ Cb {Cr — Cb ^ Ji — J2), corresponding to whether 
r' is in the red or black plaquettes of the blue lattice 
(Fig. . Retracing the steps leading to the Coulomb gas 
representation we finally obtain: 



nt^-irJ J d^r j dVn(r)n(r')ln(|r-r'|/a) 

+ E^^ J d\n^r)+E'^ J d^rn^r) , (22) 

where Ec^'''^ are the core energies of (anti)vortices located 
in red (black) plaquettes and i?^ — i?^ ~ Ji — J2. Ec^''^ 
are treated as adjustable parameters, chosen to best re- 
produce the energetics of the original Hamiltonian (|14ll . 
The explicit lattice version of H20I) follows from Ref. 1261 
where a duality transformation and a Migdal-style renor- 
malization procedure have been applied to the XY model: 

nt = ~7rjJ2s{r)s{r')H\r~r'\/a) 



i?:^.2(r)+£;,^^.2(r) . (23) 



reB 



the blue lattice in Fig. |5| Eq. (|^ can be viewed as a 
convenient lattice regularization of (|22l) . It is rederived 
in the Appendix with the help of Villain approximation. 

We can also recast Tiy in a form which interpolates be- 
tween the continuum H22() and lattice (|23|l Coulomb gases. 
For convenience, we now limit our attention to only ±1 
vortex charges (vortices and antivortices) since these are 
the relevant excitations in the pseudogap regime. The 
underlying lattice in (|23|1 is "smoothed out" into a pe- 
riodic potential V{r) whose minima are located in the 
plaquettes of the blue lattice (at vertices of the dual blue 
lattice). The precise form of V{r) can be computed in 
a continuum model of realistic cuprates involving local 
atomic orbitals, self-consistent computations of the pseu- 
dogap A, etc. and is not important for our purposes; 
only its overall symmetry matters and the fact that it is 
sufficiently "smooth" . This finally produces: 



ni = -TTj I d^r / dVn(r)n(r')ln(|r-r'|/a) 



+ / d'^rV{r)p{v) , (24) 



In H23|) s(r) = 0, ±1(±2, ±3, . . . ) and r are summed over 
the red {TV) and black {B) sites of the lattice dual to 



where n(r) is the vortex charge density as before and 
p{t) — ^ - 5{y — r") -I- Yl- 5{t — r^) is the vortex particle 
density irrespective of vortex charge. If E"^ — E^, the 
periodicity of V{y) coincides with that of the blue lattice 
- the value at the minimum is basically = E^. On the 
other hand, if E"^ ^ E^ as is the case in real cuprates, 
V{r) has a checkerboard symmetry on the blue lattice, 
with two different kinds of local extrema; the red one, at 
energy set by E^ and a black one, at energy E^ (this is 
all depicted in Fig. 

The most important consequence of these two different 
local extrema of V{r) is that there are two special loca- 
tions for the position of vortices; on the original Cu02 
lattice of Fig. [3 these correspond to vortices either re- 
siding in its plaquettes or at its vertices. On general 
grounds, we expect one of these positions to be the abso- 
lute minimum while the other assumes the role of either a 
locally stable minimum or a saddle point with an unsta- 
ble direction. Whether EJ^ < E^ or the other way around 
and whether the higher energy is a local minimum or a 
saddle point, however, can be answered with certainty 
only within a model more macroscopic than the one used 
here. In particular, a specific analysis of the electronic 
structure of a strongly correlated vortex core is necessary 
which goes far beyond the XY-type model used in this 
paper. For our purposes E^ and E^ can be treated as 
adjustable parameters and it suffices only to know that 
the microscopic physics selects either the red or the black 
plaquettes in Fig. |21 as the favored vortex core sites and 
relegates the other to either the metastable or a weakly 
unstable status. 
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IV. QUANTUM PHASE FLUCTUATIONS 

The previous discussion pertains to the thermal (or 
classical) 2D XY model. The key difference from the 
usual Coulomb gas representation for vortex-antivortex 
pairs turned out to be the inequivalence of vortex core 
positions on the blue lattice (Fig. [2Jl, which served to 
recover the translational symmetry of the original Cu02 
lattice. This effect was represented by the blue lattice 
potential ^(r) in the vortex Hamiltonian (|24() . 

To obtain the quantum version of the fluctuation prob- 
lem or an effective (2+l)D XY-type model we need to in- 
clude the imaginary time dependence in ipi — > (pi{T). In 
general, this will result in vortex-antivortex pairs prop- 
agating through imaginary (Euclidean) time, as vortex 
cores move from one plaquette to another through a se- 
quence of quantum tunneling processes. There are two 
aspects of such motion: The first is the tunneling of the 
vortex core. This is a "microscopic" process, in the sense 
that a detailed continuum model, with self-consistently 
computed core structure, is needed to describe it quan- 
titatively. For the purposes of this paper, all we need to 
know is that one of the final components characterizing 
such a tunneling process is the "mass" M of the vortex 
core as it moves through imaginary time21. This implies 
a term 

in the quantum mechanical vortex action, where r^^"''(T) 
is the position of the ith (anti)vortex in the Cu02 plane 
at time r. Note that M is here just another parameter in 
the theory which, like J, Ji, and J2, is more microscopic 
than our model^. However, the microscopic physics of 
the cuprates clearly points to M being far smaller than 
in conventional superconductors. Cuprates are strongly 
coupled systems, with a short coherence length ^ ^ kp^ 
and the vortex core of only several lattice spacings in 
size. Consequently, there are only a few electrons "inside 
a core" at any given time, making (anti)vortices "light" 
and highly quantum objects, with strong zero-point fluc- 
tuations. Furthermore, the core excitations appear to be 
gappe(& by the combination of strong coupling and local 
Mott-Hubbard correlationo^^i^°i^^i^^i^^. This has an im- 
portant implication for the second component of the core 
tunneling process, the familiar Bardeen-Stephen form of 
dissipation, which is such a ubiquitous and dominating 
effect in conventional superconductors. In conventional 
systems the core is hundreds of nanometers in diameter 
and contains thousands of electrons. M is large in such 
a superconductor and, as a vortex attempts to tunnel to 
a different site, its motion is damped by these thousands 
of effectively normal electrons, resulting in significant 
Bardeen-Stephen dissipation and high viscosity. The mo- 
tion of such a huge, strongly damped object is effectively 
classical. Vortices in under doped cuprates are precisely 
the opposite, with small M , the Bardeen-Stephen dissi- 



pation nearly absent and a very low viscosity. In this 
sense, the quantum motion of a vortex core in under- 
doped cuprates is closer to superfluid ■*He than to other, 
conventional superconductors^^. We are therefore justi- 
fied in assuming that it is adequately described by H25() 
and ignoring small vortex core viscosity for the rest of 
this paper—. The special feature of the cuprates, how- 
ever, is the presence of gapless nodal quasiparticles away 
from the cores which generate their own peculiar brand of 
(weak) dissipation - such effect is discussed and included 
later in this section. 

The second aspect of the vortex propagation through 
imaginary time involves the motion of the superflow ve- 
locity field surrounding the vortex outside its core, i.e., in 
the region of space where the magnitude of A is approx- 
imately uniform. This is a long range effect and, unlike 
the vortex core energies, mass, or Bardeen-Stephen dissi- 
pation, exhibits certain universal features, shared by all 
superconductors and superfluids. For example, in super- 
fluid "^He this effect would result in a Magnus force acting 
on a vortex. In a superconductor, the effect arises from 
the time derivative of the phase in the regions of uniform 
A. The origin of such time derivative and its physical 
consequences are most easily appreciated by considering 
an ordinary (fictitious) strongly fluctuating s-wave super- 
conductor with a phase factor exp(z0i(r)) at each site of 
the Cu02 lattice. For reader's benefit we discuss this 
case first as a pedagogical warm-up for what lies ahead. 



A. Pedagogical exercise: quantum fluctuating 
s-wave superconductor 

After performing the FT transformation and forgetting 
the double valuedness problem (ignoring the Berry gauge 
field a) since we are concentrating on the charge channel, 
the quantum action will contain a purely imaginary term 
(see Section II): 

-/ dTV(n,;t(T)-fn,^(T))9,0, , (26) 
^ Jo i 

where riia- = ipia-ipicr. It is useful to split the electron 
density into its average and the fluctuating parts: n = 
n^' + ni n + 6n. The average part is unimportant 
for the spin- wave phase due to the periodicity of ^^(t) 
in the interval [0,/3]. In the vortex part of the phase, 
however, this average density acts as a magnetic flux 
seen by (anti)vorticesi^. This first time derivative is just 
the charge Berry phase and it couples to (half of) the 
total electronic charge density. After the fermions are 
integrated out, the fluctuating part of the density will 
simply generate a 4>f term whose stiffness is set by the 
fermionic compressibility. 

For simplicity, we will first ignore the charge Berry 
phase and set h — (mod 2). This results in a (2-1- 1)D 
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XY model: 

/:jfy = ^A'o^</'- - J^cos(0.-<^,) , (27) 

i nn 

where Kq is effectively the fermionic charge susceptibility 
(compressibility): Kq ~ Xc ~ i^n-'i)- Again, we expand 
the cosine to the leading order (the same results are ob- 
tained in the Villain approximation, see Appendix) and 
separate the regular and singular (vortex) contributions 
to 9^0 — > dfj,x + {dfj.v)v As a result, Cxy is transformed 
into 

^K^id^cj,)' ^ iW^d,x + ^W^id^cly), + ^K;;'W^ , (28) 

where Kf^ = {Ko,K,K) and is a real Hubbard- 
Stratonovich three-vector field. The integral over the 
free field x can be carried out, producing the constraint 
d ■ W = 0. The constraint is solved by demanding that 
W = d X Ad, where will soon assume the meaning of 
the dual gauge field. The vortex part is now manipulated 
into: 

Cxy -> i{d x A^) ■ (90). + Ik'^O x Ad)l 

^-2mAd-n+^K;^HdxAd)l , (29) 

where the integration by parts and d x {d(t>)y = 27rn have 
been used and is the vorticity in the (2-1-1) dimen- 
sional spacetime (vorticity, by its very nature, is con- 
served: d ■ n = 0). Now, we use the standard transi- 
tion from Feynman path integrals to coherent functional 
integrals^. A relativistic vortex boson complex field 
^(x) is introduced, whose worldlines in (2-1-1) dimen- 
sional spacetime coincide with fluctuating vortex loops 
(see the Appendix for a more detailed account). The 
first term in 1)29(1 is just the current of these relativistic 
vortex particles coupled to a vector potential Ad- Fur- 
thermore, vortices have some intrinsic line action J dsSo, 
coming from core terms (and/or lattice regularization) 
which supply the (2-1-1) kinetic term. In the end, one 
obtains a dual Lagrangian: 

£d = \id+t2nAdm'+m'{rm' + ^m^+^^idxAd)l 

(30) 

where g > is a short range repulsion describing the 
penalty for vortex core overlap. The mass term will in 
general be spatially modulated, reflecting the underlying 
lattice potential, as in (|24|l . A detailed derivation of the 
above dual Lagrangian can be found in the Appendix. 

It is useful to underscore the following relation between 
the steps that led to Ij^UI) and the formalism discussed in 
Section II, which is coached in the language closer to 
the original electrons: Only the familiar long range Biot- 
Savart (anti)vortex interactions, mediated by dual gauge 
field ^rf, arise from the fermionic action in Section II and 
Refs. UalsM This corresponds to large regions of space 



where the pseudogap A is large and approximately uni- 
form. In contrast, small (anti)vortex core regions, where 
A might exhibit significant non-uniformity, supply the 
the core kinetic part (|25|l . vortex core energy and the 
short range core repulsion term (g), which are all stored 
in the "Jacobian" for v and a gauge fields. This is signif- 
icant since it enables us to deal rather straightforwardly 
with the charge Berry phase which now must be restored. 

First, notice that the self-action of Ad, which was intro- 
duced in this section through a Hubbard-Stratonovich de- 
coupling in Eq. (|28|l . actually follows "microscopically" 
from the integration over the fermions and the gauge field 
V in Section II andpi^. There Ad was introduced as a field 
enforcing a 5-function constraint i5[(9 x ■i;)/7r — n] but the 
physics is precisely the same, the formal difference being 
just the order in which we integrate the fermionic matter. 
Once we restore the Berry phase term (i/2) / (fixn{dr4>)v 
back to the action, we note that in the formalism of Sec- 
tion II and^- none of the fields in the dual Lagrangian 
(|Sn|l couples to it directly. Instead, it is vq that enters the 
Berry phase via the i5-f unction constraint dr4> vq. The 
Berry phase term will then affect (|30|l via the coupling 
of vo to the spatial part of Ad, i.e., the dual magnetic 
induction = V x A^. 

The above observations suggest that it useful to sepa- 
rate out the saddle point part of w as: v ~f ~iv + Sv. 
The saddle point part v is determined by minimizing 
the total action (the "— " sign is chosen so that vq cou- 
ples to the electron density as a chemical potential). Of 
course, if there were no Berry phase term we would have 
-0 = 0. With the Berry phase term included v = but 
Vq is now finite. Observe from Eqs. (1,2) of Ref. 
that the saddle point equation for vq simply reduces to: 
in = Brf = V X Ad (note that, relative to our notation, 
dual gauge field in Ref. is rescaled as Ad 2'jTAd)- 
So, with the Berry phase included, the final form of the 
dual theory still remains given by Eq. ((30|l but must be 
appended by the constraint ^{n{x)) — (Bdix))- 

If we now apply the dual mean-field approximation to 
(|28|l one obtains J-mf given by: 

|(V + z2^A,)<i>p + m2(r)|$|2 + + ,(31) 

with the constraint taking the form Brf(r) = i(n(r)). 
One immediately recognizes (|31|l as the Abrikosov- 
Hofstadter problem for a dual type-II superconductor 
- we are in type-II regime^ since small compressibility 
implies large dual Ginzburg parameter Kd ^ ^/VKo ^ 
1/y/Xc ~ subjected to the overall (dual) magnetic field 
flux per plaquette of the dual lattice, /, given by / = 
p/q ~ (1 — x)/2, where x is doping. Note that at half- 
filling h = 1 ^ f — 1/2. The solution that minimizes 
(|31|l at half-filling is a checkerboard array of vortices 
and antivortices in with the associated modu lation 
in B(j(r) = i(n(r)), as discussed in Refs. Il2l38i This 
is nothing but the dual equivalent of a CPCDW in an 
s-wave superconductor at half-filling, with a two-fold de- 
generate array of alternating enhanced and suppressed 
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charge densities on sites of the original lattice. An over- 
whelming analytic and numerical support exists for this 
state being the actual ground state of the negative U 
Hubbard model^^, the prototypical theoretical toy model 
in this category. This gives us confidence that the mean- 
filed approximation captures basic features of the prob- 
lem. 



B. Quantum fluctuating d-wave superconductor 

The above pedagogical exercise applies to an idealized 
strongly fluctuating s-wave superconductor and its asso- 
ciated ordinary (2-1- 1)D XY model. What about our d- 
wave lattice superconductor (LdSC) and its effective XY- 
type model H14(l . with the bond phases mapped into site 
phases on the blue lattice? Again, we go from exp{i9ij) 
to exp{i(pk) as before. The cosine part of H27(l is replaced 
by TCxY 114|l and is handled in the same way as for the 
classical case. This will ultimately result in a modulated 
potential V{r) of Eq. (|24l) . In the context of the dual 
theory H30() this will translate into a position-dependent 
mass term m^(r)|$(a;)p in the dual Lagrangian H3U|) ex- 
cept now this mass term has the checkerboard symmetry 
on the blue lattice, with rn^ in the red plaquettes differ- 
ent from ml in the black plaquettes {nri^ — ml ^ E'^ — E^). 
Including the constraint on the overall dual induction 
seems to complete the process. 

Alas, the situation is a bit more involved. First, the 
time derivative part and the Berry phase are more com- 
plicated. The difficulty arises since the bond supercon- 
ducting phase 0ij of a LdSC couples in a more com- 
plicated way to the site fermionic variables. However, 
we can still deal with this by enlisting the help of the 
FT gauge transformation which screens the long-range 
part of 6ij by + l^j: where (f>i are the phases in the 
fermionic fields. After this transformation (we are again 
ignoring which we can put back in at the end) one 
obtains i\nidr(f)i at each site of the Cu02 lattice. This 
translates into ^nidr(l>i + '^n'j^r't'j for each bond (ij) of 
the Cu02 lattice. This bond expression can be rewritten 
as: 

i 1 i 1 

-[rii + nj]-{dr4>t + dr4>j) + ^j]2^(^r(f>i - dr4>j) 

= + + (•••) , (32) 

8 

where (■ • • ) contains higher order derivatives and is typi- 
cally not important in the discussion of critical phenom- 
ena (but see below). In the end, following our replace- 
ment of bond phases on Cu02 lattice with site phases 
on the blue lattice Oij — > (/Sfc we finally obtain the Berry 
phase term of our lattice d-wave superconductor: 

z X % 

-[n^ + nj]dripk o + "jl^Tiy^fc + -SAijdrifk ■ (33) 

o o o 

The expression ()32f) seems too good to be true and it 
is - the replacement 9ij + holds only far away 



from vortices, when the phases change slowly between 
nearby bonds (sites). In general, the bond phase of the 
LdSC will not couple to the electron densities in a simple 
way suggested by (I32|l . Still, we have gained an impor- 
tant insight; its coupling to the overall electron density 
represented by the leading term in H33|l is effectively ex- 
act. This follows directly from the electrodynamic gauge 
invariance which mandates that both regular and singu- 
lar configurations of the phase must have the same first 
time derivative in the action. For arbitrarily smoothly 
varying phase, with no vortices present, the replacement 
dij ^i'i ~^ \4>3 is accurate to any desired degree and 
the Berry phase given by the leading term in H33|l follows. 
Incidentally, this result is not spoiled by the higher or- 
der terms, represented by (• • • ) in (|32(l . since they do not 
contribute to the overall dual magnetic field, by virtue of 
fii — fij = 0. We, however, cannot claim with the same 
degree of certainty that (5Aij is in a simple relationship 
to variations in the electronic densities on sites i and j 
as would have ili^. Instead, we should view 5/S.ij 
in (|^ as a Hubbard-Stratonovich field used to decouple 
the 9^^ term in the quantum action. This makes 5l\ij in- 
herently a bond field whose spatial modulation translates 
into variation in the pseudogap A and will be naturally 
related to the dual induction below^. This is dif- 
ferent from our pedagogical s-wave exercise where the 
modulation in B^ was directly related to the variation in 
the electronic density. 

The second source of complication is more intricate 
and relates to the fact that, as T ^ 0, we must be con- 
cerned about nodal fermions. In the s-wave case and the 
finite temperature d-wave XY-type model (|14|) we know 
that coupling constants J, Ji, J2, and so on, are finite 
or at least can be made so in simple, reasonably realistic 
models. In general, we need to keep only several of the 
nearby neighbor couplings to capture the basic physics. 
There is a tacit assumption that such expansion is in fact 
well-behaved. In the quantum d-wave case, however, the 
expansion in near neighbor couplings is not possible - 
at T = its presumed analytic structure is obliterated 
by gapless nodal fermions. Naturally, after the electrons 
are integrated out we will still be left with some effective 
action for the phase degrees of freedom 9ij (r) but this ac- 
tion will be both non-local and non-analytic in terms of 
differences of phases on various bonds. The only option 
appears to be to keep the gapless fermionic excitations 
in the theory. This in itself of course is perfectly fine 
but it does not advance our stated goal of "bosonizing" 
the CPCDW problem a least bit. Note that the issue 
here is not the Berry gauge field a and its coupling to 
nodal fermions - even when we ignore coupling of vor- 
tices to spin by setting a ^ and concentrate solely 
on the charge sector as we have done so far, the gapless 
BdG quasiparticles still produce non-analytic contribu- 
tions to the phase-only action. In short, what one has 
here is more than a problem; it is a calamity. The im- 
plication is nothing less but that there is no usable XY- 
type model for the charge sector of a quantum d-wave 
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superconductor, the optimistic title of the present paper 
notwithstanding. 

Fortunately, there is a way out of this conundrum. 
While indeed we cannot write down a simple XY-like 
model for the fluctuating phases it turns out that the 
dual version of the theory can be modified in a relatively 
straightforward way to incorporate the non-analyticity 
caused by nodal fermions. To see how this is accom- 
plished we go back to Section If and use our division 
of electronic degrees of freedom into low-energy nodal 
{ipcr,a) and "high-energy" anti-nodal fermions, tightly 
bound into spin-singlet Cooper pairs {'4>a,{ai3))- Imagine 
for a moment that we ignore nodal fermions altogether, 
by, for example, setting the number of Dirac nodal flavors 
A'' ^ (in a single Cu02 plane N — 2). Exact integra- 
tion over anti-nodal electrons i^^p) leads to an XY-type 
model precisely of the type discussed at the beginning 
of this subsection. Since anti-nodal electrons are fully 
gapped the effective Hamiltonian has the conjectured 
form p4|l with finite couplings J, Ji and J2, and with 
higher order terms which decay sufficiently rapidly. The 
explicit values of J, Ji and J2 at T = might be difficult 
to determine but not more so than for an s-wave super- 
conductor. Furthermore, since anti-nodal fermions carry 
all the charge density, the Berry phase term remains as 
determined earlier and so does the timelike stiffness of 
ipi. Derivation of the dual theory proceeds as envisioned, 
with all the short range effects stored in vortex core terms 
and ultimately in m^(r), and with only remaining long 
range interactions mediated by the dual gauge field Ad- 
We now restore nodal fermions {N 7^ 0); as already em- 
phasized, the only coupling of nodal fermions to the dual 
theory is via the gauge fields v and a. Since we are ig- 
noring the spin channel we can neglect a. As far as v is 
concerned, including nodal fermions leads to a non-local, 
non-analytic correction to its stiffness, expressed as 

where is determined by anti-nodal electrons. Qfj,u{q) 
is the contribution from nodal fermions, linear in q 
and in general quite complicated. We give here few 
simpler limits: (5oo(0,w„) = c|w„|, Qo* = QiO = 0, 
Qxx{ci,0) = Qyy{(l,0) = -cg>, Qxy{q,0) = Qyx{ci,0) = 
-cq<sgn-{qx,qy), where g> = maxdg^;], g< = 

min(|g^|, \qy\) and c - ttVIGv^. 

In the language of dual theory, this translates into a 
modified self-action for the dual gauge field: 

-^^{dx Ad)^[K^S^, + NQ^,]-\dx Ad). , (34) 

where [Kf^Sfj,. + NQf^.]^^ is the tensor inverse of Kfj,S^,j + 
NQfiu- This action is clearly more complicated than its 
Maxwellian s-wave counterpart but nevertheless decid- 
edly manageable - most importantly, the part of the 
dual action involving dual boson field $ remains unaf- 
fected. Note that nodal fermions induce subleading but 



still long-ranged interactions between vortices, in addi- 
tion to the familiar Biot-Savart interactions. These in- 
teractions have square symmetry on the blue lattice, re- 
flecting their nodal origin. Furthermore, when Ad is in- 
tegrated over, they will produce a peculiar dissipative 
term ^ \ujn\vo which describes the damping of collective 
quantum vorticity fluctuations by gapless nodal quasipar- 
ticles. The importance of this effect is secondary relative 
to the mass term for vq which is always generated by anti- 
nodal fermions, basically because the density of states 
N{E) for nodal fermions vanishes as ^ 0. However, 
given the smallness of traditional Bardeen-Stephen core 
viscosity this "nodal" mechanism is an important source 
of dissipation of vortex motion in underdoped cuprates. 
Since our focus in the present paper is dual mean-field 
theory, this dissipative term will play no direct role. 

Armed with the above analysis we can finally write 
down the Lagrangian of the quantum XY-type model de- 
scribing a LdSC: 

^XY — * ^ fi'Pi + 5^ V'l + '^XY + ^nodal + ^coro 

i i 

i i nn 

- ^ C0s((^i -ipj) - J2 ^ COs{ipi - (fj) 
rnnn bnnn 

+ >C„odal[cOs((/?j - + £corc , (35) 

where the sums run over sites i of the blue lattice and 
the notation is the same as below Eq. (|14|) . fi — 
I (fife -I- fij) with k,j being the end sites of bond i, the 
time-like phase stiffness Kq results from the Hubbard- 
Stratonovich integration over SAij as detailed above, and 
'Ccoic contains core terms coming from a small region 
around the (anti)vortex where A itself is significantly 
suppressed, an example being 1)25(1 . The explicit form 
of £nodai[cos(c/Ji — ipj)] in thc XY model language is un- 
known but it is a non-analytic, non-local functional of 
cos{ipi — <y3j); its effect will be incorporated once we arrive 
at the dual description, following the above arguments. 
Note that the short range parts of cosine functions and 
time derivative in (|35|l will be subsumed into £core once 
we go to continuum or Villain representations of the prob- 
lem, as described elsewhere in the paper. Furthermore, 
observe that at half-filling the dual flux per plaquette of 
the blue lattice has an intrinsic checkerboard pattern: 
/= 1/2 for black plaquettes and / = for red plaque- 
ttes. This is a direct consequence of (fir) being a bond 
phase, residing on the blue lattice sites in Fig. [3 its Berry 
phase given by Eq. H33|l . Thus, in the quantum problem, 
the Berry phase combines with the next-nearest neighbor 
bonds (for Ji ^ J2) in breaking translational symmetry 
of the blue lattice down to the checkerboard pattern of 
Fig. 12 The average flux through the unit cell on the 
blue lattice is / = 1/2 which is just as it should be since 
there are two plaquettes of the blue lattice per single pla- 
quette of the original black lattice and the dual flux per 
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plaquette of the Cu02 lattice is / = 1/2 at half-filling. 

It is now straightforward to derive a dual representa- 
tion of H35(l by retracing the path that led from (|27|l to 
(|3U|I . In this fashion we obtain the dual Lagrangian of the 
(2+l)D XY model appropriate for a strongly fluctuating 
d-wave superconductor near half-filling^: 

+ \{d X Ad),\K^6,,, + NQ^,]-\d X Aa), + {■■■) , 

(36) 

where d x Ad satisfies the constraint arising from the 
charge Berry phase. The dual Lagrangian for a d-wave 
superconductor (|36|l differs from its s-wave counterpart 
(|30|l in the following important respects: i) the period- 
icity of the mass potential is different and reflects the 
checkerboard periodicity on the (dual) blue lattice set 
by our vortex core potential V{y), ii) the self-action for 
the dual gauge field contains non-local, non-analytic con- 
tribution of nodal fermions, Hi) the modulation of the 
dual flux Bd relates to the variation in the pseudogap 
5Aij (as opposed to the electronic density 5ni) via the 
constraint ^B(i(r) ^ (l/16)(5Aij(r)), and iv) at dop- 
ing x, the overall dual flux through the unit cell of the 
blue lattice imposed by the Berry phase constraint is 
/ = p/q = (1 — x)/2, comprised of / = p/q = (1 — x)/2 
at a black plaquette and / = at a red one. We there- 
fore expect the results for stable CPCDW states to be 
different from those of a hypothetical strongly fluctuat- 
ing s-wave superconductor. Also, note that Kq is now 
set by Kq ^ (jJAyp), which is still related to fermion 
compressibilit}*^ and still relatively small since the ba- 
sic aspect of our approach is that amplitude fluctuations 
in Aij are suppressed. This implies our dual supercon- 
ductor remains in the type-II regime {na ~ ^1). 
Finally, (• • • ) denotes higher order terms which have been 
neglected. For those readers who find the above road to 
perhaps a bit too slick, we give a detailed step-by- 
step derivation in the Appendix. 

An important point should be noted here: The ne- 
glected higher other terms (■ • • ) involve higher order ki- 
netic terms, additional intermediate range interactions, 
short range current-current interactions and numerous 
other contributions. All are irrelevant in the sense of 
critical behavior. In our problem, however, we are not 
interested only in the critical behavior. In particular, 
we would like to determine the charge modulation of the 
CPCDW on lengthscales which are not extremely long 
compared to the Cu02 lattice constant. These additional 
terms could be important for this purpose. Particularly 
significant in this respect arc the higher order corrections 
to the Berry phase (|32|) which, while not affecting the 
value of overall /, do influence the form of the effective 
dual Abrikosov-Hofstadter problem that ensues. 

The strategy of keeping a large number of otherwise ir- 
relevant terms is not a practical one. We will therefore in- 
troduce a simplification here which is actually quite nat- 



ural and allows us to retain the essential realistic features 
of the original model and maintain the particle-hole sym- 
metry as well. Consider a situation where m^(r) in H3t)|) is 
very strongly modulated. This is a "tight-binding" limit, 
universally considered appropriate for cuprates, and we 
can simply view quantum (anti)vortices as being able to 
hop only between the nearby plaquettes of the blue lat- 
tice, as is clearly implicit in (|35|l . Furthermore, we as- 
sume that the two extrema in (r) are separated by an 
energy scale larger than such hopping. This is a perfectly 
reasonable assumption since it illustrates an already im- 
portant characteristic of our effective model, the fact that 
red and black plaquettes of the blue lattice are intrinsi- 
cally not equivalent (by the form of the Berry phase and 
Ji 7^ J2). Under these circumstances we can rewrite the 
dual Lagrangian in a tight-binding form: 



= ^ 1(9. + i2TiAM)'^.r? + ^ 1(9. + i2nAdM'' 

r b 




r b 



+ -{d X Ad)^[K^Sf,, + NQ^,]-\d X Ad),. + (•••) , 

(37) 



where vortex boson fields ^r(b) reside on red (black) 
plaquettes of the blue lattice, irr(&) is the vortex 
hopping between the nearest red-red (black) neighbors, 

exp(i27r ^''^ ds ■ Ad) are the appropriate Peierls factors, 
— TO^I ^ trb (core energies on red and black plaque- 
ttes are significantly different), and V x in the last 
term is the lattice flux of the dual magnetic field, equal 
/ = p/q = (l~2;)/2 per unit cell of the original Cu02 lat- 
tice. We have also assumed that it is the red plaquettes 
that are favored by vortex cores, making it unnecessary 
to include ttb explicitly. The assumption Er <C E(, ap- 
pears naturally warranted by the overall symmetry of the 
problem but, should the details of microscopic physics 
intervene and reverse the situation in favor of the black 
plaquettes, all one needs to do is exchange labels r 6 
(note that < 0, rri^ > in a dual superconductor). 
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The resulting tight-binding dual Hamiltonian: 
Hd = - ^ Ur exp{i2TT I ds-Ad)$*$r' 

(rr') 

- ^trbex.p{i2Tr ds • Ad)^>*$6 - (c.c.) + 

(rb) 
r 

+ + + (...) (38) 

b 

will be analyzed in two limits: for trr ^ i^b/l^r ~ 
in which case we can simply set trb — + 0, the effective 
Hofstadter problem assumes the form equivalent to the 
standard s-wave case defined on the sites of the red lat- 
tice (Fig. I2J). Importantly, however, the relation between 
the modulation in and SAij remains different from 
the s-wave case and peculiar to a d-wave, as explained 
above. We will call this limit an HI model. Similarly, 
in the opposite case ^ ^r&ZI'^r ~ ™&l can set 
trr ^ and obtain an effective Hofstadter problem with 
hoppings between red and black sites only. Note that this 
situation is not equivalent to the standard s-wave case: 
In order to hop from one red plaquette to another a vor- 
tex must go through a black site, picking up a Peierls 
phase factor different from the one for a conventional 
direct red-to-red hop. In considering these Peierls fac- 
tors we must exercise caution since the "assisted" hops 
between red plaquettes pass directly through the dual 
fluxes } — pI q — (1 — x)/2 located on black sites. By 
infinitesimally displacing the said flux, one is back to the 
situation where all closed paths of hops are composed of 
fluxes / = p/q = (1 — x)/2 through black and / = 
through red plaquettes of the blue lattice. The resulting 
Hamiltonian has an exact particle-hole symmetry around 
half-flUing (x —x), as it should. We call this the H2 
model. Finally, still in the limit trr ^ irb/l™r~™6l with 
we can "spread" the dual flux so it is uniformly 
distributed throughout each blue lattice plaquette and 
equal / = p/q = (1 — a;)/4. In this situation, dubbed an 
H3 model, the x — > —x symmetry is obeyed only approx- 
imately, for small x (near half-filling), but this is all we 
are interested in. We have established by explicit com- 
putations that the H3 model satisfies the particle-hole 
symmetry at some Hofstadter fractions / while it does 
not appear to do so at others; it is for this reason that 
we tend stay away from the H3 model in this paper. The 
reader should note that the issue of how to deal with dual 
fiuxes when hopping through black plaquettes arises only 
in tight-binding models H2 and H3 (but not HI!) - the 
original dual Lagrangian (|36|) is free of such issues and has 
manifest x — > —x symmetry. The down side, of course, 
is that such a continuum theory is far more difficult to 
solve, both analytically and numerically. 

The Hamiltonian (|38|l . in its three editions HI, H2, and 
(to a limited extent!) H3, defines probably the simplest 



dual version of the Hofstadter problem which appropri- 
ately builds in the d-wave character of the fiuctuating 
lattice superconductor and the essential phenomenology 
of the pseudogap state in underdoped cuprates. The ex- 
tra terms (• • • ) can still be used for fine-tuning (for ex- 
ample, g gr,gb, additional |$rP|*i'bP repulsions, etc.) 
but the important particle-hole symmetry around half- 
filling is already present without them. All the detailed 
numerical calculations reported in the latter pages of this 
paper and described in Ref. use the mean-field ver- 
sion of and the three models based on it, HI {trr ^ 

tlb/l'^r-^bl) and H2 and H3 {trr < irb/l"^r ^ "^fcl)' as 
the point of departure. 



V. DUAL SUPERCONDUCTOR AND COOPER 
PAIR CDW IN UNDERDOPED CUPRATES 

The previous sections concentrated on the derivation 
of the effective quantum XY-type model for phase fluctu- 
ations in underdoped cuprates and its dual counterpart 
and (|S7|) . In this section we pause to take stock of 
where we are with respect to the real world and to con- 
sider some general features of the physical picture that 
emerges from Eqs. (I3m37|) . First, the dual supercon- 
ductor H36() describes the physics of strongly fluctuating 
superconductors in terms of a complex bosonic fleld $, 
which creates and annihilates quantum vortex-antivortex 
pairs viewed as "charged" relativistic particles (vortices) 
and antiparticles (antivortices) - this is depicted in Fig. 
O The conserved "charge" is just the vorticity associated 
with these topological defects, -1-1 for vortices and -1 for 
antivortices, and the gauge field coupled to it is nothing 
but Ad- Physically, Ad describes the familiar logarithmic 
interactions between (anti)vortices. We stress that the 
dual description is just a convenient mathematical tool: 
Its main advantage is that it allows a theorist to access 
a strongly quantum fluctuating regime of a superconduc- 
tor, where the superfluid density ps at T = can be very 
small or even vanish, while the pairing pseudogap A re- 
mains large. According to the theory of Ref. this is 
the regime that governs the properties of the pseudogap 
state in underdoped cuprates. This regime is entirely 
inaccessible by more conventional theoretical approaches 
which use the mean-field BCS state as the starting point 
around which to compute fluctuation corrections. 

The dual superconductor description predicts two 
basic phases of cuprates: ($) = is just the fa- 
miliar superconducting state. Quantum and thermal 
vortex-antivortex pair fluctuations are present (and thus 
(|$|2) ^ 0) but these pairs are always bound, result- 
ing in a flnite, but considerably suppressed superfluid 
density. As vortex-antivortex pairs unbind at some dop- 
ing X — Xc, the superfluid density goes to zero and the 
superconducting state is replaced by its dual counter- 
part, ($) 7^ 0. The meaning of dual ODLRO is actually 
quite simple: Finite ($) means that vortex loops (loops of 
vortex-antivortex pairs being created and annihilated in 
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FIG. 3; Quantum fluctuations of vortex-antivortex pairs, xy- 
plane is the Cu02 layer and r-axis shows direction of imag- 
inary time. Vortices (arrows pointing upwards) and antivor- 
tices (arrows pointing downwards) are always created and an- 
nihilated in pairs. Note that the structures arising from linked 
creation-annihilation events form oriented loops carrying ±1 
vorticity. These loops are just the virtual particle-antiparticle 
creation and annihilation processes in the quantum vacuum of 
the relativistic bosonic field $(a;), as described in the text and 
the Appendix. $ is our duaZ order parameter: In a physical su- 
perconductor such vorticity loops are finite and ('!>) = 0. Note 
that the intersections of such finite loops with the xy-plane 
at any given time r define a set of bound vortex-antivortex 
dipoles in the Cu02 layer. The superconducting order is lost 
when vortex-antivortex pairs unbind and the average size of 
the above vorticity loops diverges - some of the loops become 
as large as the system size. This is the pseudogap state, with 
($) 7^ 0. The reader should note the following amusing aside: 
The above figure can easily be adapted to depict the low en- 
ergy fermionic excitations of theory The creation and 
annihilation processes now describe spin up (arrows pointing 
upwards) and spin down (arrows pointing downwards) quasi- 
particle excitations from the BCS-type spin-singlet vacuum. 
The loops carry a well-defined spin and can be though of as 
relativistic BdG Dirac particles/antiparticles - their massless 
character in a nodal d-wave superconductor implies the pres- 
ence of loops of arbitrarily large size. These fermionic loops 
move in the background of fluctuating vortex loops discussed 
above, their mutual interactions encoded in gauge fields v and 
a. This is a pictorial representation of the theory iQ. 



(2-t-l) dimensional spacetime) can now extended over the 
whole sample, i.e. the worldline of a dual relativistic bo- 
son can make its way from any point to infinity (see Fig. 
IHl . The presence of such unbound vortex-antivortex ex- 
citations directly implies vanishing helicity modulus and 
thus the absence of the Meissner effect and ps = (see 
Section II for details). The phase diagram of a dual su- 
perconductor as it applies to cuprates is shown in Fig. 

m 

In the dual mean-field approximation, just as in a con- 
ventional one, we ignore fluctuations in $ and minimize 
the action specified by H36|l with respect to a complex 
function $(r). In a physical superconductor {x > Xc) 
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FIG. 4: The schematic pha se d iagram of underdoped 
cuprates in the theory of Refs. I5ll2l T* denotes the pseu- 
dogap temperature (T* ~ A). Dual superconducting order 
(finite {$)) implies the absence of the true, physical super- 
conductivity. Conversely, the absence of dual order ({$) — 0) 
corresponds to the physical superconductor. The shaded area 
represents the region of coexistence between a strongly fiuc- 
tuating superconductor and a CPCDW state, which will gen- 
erally occur in the dual theory. This region of "supersolid" 
behavior is characterized by {"!>) — but finite (I'l'l^), which 
is modulated within the Cu02 layer according to our theory 
of a dual superconductor. The precise size of a coexistence 
region, however, is difficult to estimate from the mean-field 
theory and a more elaborate approach, including fiuctuations 
in $ and Ad, needs to be employed. Finally, the finite tem- 
perature phase boundary of the CPCDW (dashed-dotted line) 
should not be taken quantitatively apart from the fact that it 
is located below T*. 



^(r) = 0. For x < Xc, the mininmm action corresponds 
to finite <I>(r). However, since the charge Berry phase 
translates into a dual magnetic field this finite <I>(r) 
must contain Nd vortices, where Nd is the number of 
the dual flux quanta piercing the system. Note that 
this number is nothing but half of the total number of 
electrons Nd — N/2, the factor of one half being due 
to the fact that we are considering Cooper pairs. In a 
tight binding representation l|37|l this implies a dual flux 
f ^ p/q = (1 ~ 2;)/2 per each Cu02 unit cell. The pres- 
ence of such vortex array in $(r) will be accompanied 
by spatial variation in Bd{r), which translates into the 
modulation of the pseudogap 6Aij and thus into a CDW 
of Cooper pairs (CPCDW). Consequently, in the vocab- 
ulary of the dual mean-field approximation, the question 
of formation and the structure of the CPCDW is equiv- 
alent to finding the Abrikosov vortex array on a tight 
binding lattice, i.e. the Abrikosov-Hofstadter problem 
defined by (|37I38|) . The formation of CPCDW results in 
a modulation of the local tunneling DOS and one is led 
to identify the "electron crystal" state observed in STM 
experiment3Si2iia as the CPCDW. 
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Our task is to determine the specific pattern in which 
the vortices in <&(r) arrange themselves to minimize the 
expectation value of the dual Hamiltonian 1)38(1 . Once 
this is known we can determine the modulation in the 
dual induction and the structure of the CPCDW fol- 
lows from the Bd ^ 6Aij correspondence. In this sec- 
tion we are interested in general qualitative results and 
we therefore focus on the HI model where such results 
are more transparent. The H2 and H3 models turn out 
to be more opaque and have to be studied numerically 
as soon as one is away from half- filling. 

The solution of the Abrikosov-Hofstadter problem is 
obtained as follows: one first sets g = in (|38() and finds 
the ground state of the resulting quadratic Hofstadter 
Hamiltonian Tidig = 0) given by 

{rr') (rb) 

-(c.c.) + J2ml\'fr\'+Y.ml\^t\' ■ (39) 

r b 

At flux f = p/q the ground state is g-fold degenerate 
and we denote it by $^*)(r). One then turns on finite 
g, forms a linear combination of these degenerate states 
oiq^'^'^'' (r) , and determines variationally the set of co- 
efficients {aq} which minimizes the full Hamiltonian (|38|l . 
With {aq} fixed in this fashion the only remaining degen- 
eracy in the ground state consists of lattice translations 
and rotations. Once the ground state is found, the 
magnetic induction B^ follows from Maxwell equation: 

^^^Ax^B,(r)-j*=0 , (40) 
oAd Kq 

where £d is given by Eq. 1(3 7() . STid = A x A, and j* is 
the current in the ground state of dual Hamiltonian l|39|) 
with the uniform dual flux /. All quantities in H4U|I are 
defined on the blue lattice of Fig. |21 A is a lattice deriva- 
tive, j* and Ad are link variables, and B^ = A x A^ is a 
site variable. The detailed definitions of all these objects 
are given in the next section. The non-local, non-analytic 
self-energy for the dual gauge field in l(T7|) was replaced 
by an effective Maxwellian {Kf^ K/^i) - this approx- 
imation is valid for weak modulation JB^. We should 
stress that this way of determining the ground state 
and dual induction Brf(r) is valid only if Kq is sufficiently 

small so that the dual Ginzburg parameter 1 / V Ko 
is sufficiently large - in underdoped cuprates this is a 
justified assumption since strong Mott-Hubbard correla- 
tions strongly suppress all charge density fluctuations. 
At higher modulation, i.e. for intermediate values of 
Kd, the nodal contribution becomes more significant and 
its intrinsic square symmetry (see discussion around Eq. 
(|34|l ') will act to orient 5'Bd relative to the blue lattice. 
If this is the case, the interplay between this "nodal" ef- 
fect and the one arising from the Abrikosov-Hofstadter 
problem itself can lead to interesting new patterns for 
the CPCDW state; a detailed study of such an interplay 



is left for the future. Finally, from (5Bd(r) = B(;(r) — / 
we obtain (5Ay which can be fed back into the electronic 
structure via the expressions given in Section II. 

Notice that the above method of solving the problem 
corresponds to the strongly type-II regime of a dual su- 
perconductor {Kd » 1). In this Hmit, the CPCDW pat- 
tern is primarily given by the dual supercurrent j* in 
the Maxwell equation 1(401) . itself determined by the so- 
lution to the Abrikosov-Hofstadter problem l(39() . The 
modulation in 5'Bd only reflects this pattern of vortices 
in j* ((40(1 . and the dual magnetic energy is only a small 
fraction of the Abrikosov-Hofstadter condensation en- 
ergy. Imagine now that we ask the following question: 
What are the interactions among vortices in $ that have 
resulted in <I''^*'^ being the ground state of the Abrikosov- 
Hofstadter problem? This question is analogous to the 
one inquiring about the interactions that have led to 
the triangular lattice of vortices in the continuum ver- 
sion of the problem, i.e. the interactions inherent in 
the famed Abrikosov participation ratio (3a- In both 
cases, these are far from pairwise and short-ranged - they 
are in fact intrinsically multi-body interactions, involving 
two-, three-, and multiple-body terms, all of compara- 
ble sizes and all long-ranged^. They can be thought 
of as the interactions among the center-of-mass coordi- 
nates of Cooper pairs. This should be contrasted with 
the picture of real-space pairs, interacting with some sim- 
ple, pairwise and short-ranged interactions. The pair 
density-wave patterns in this case are determined not 
by the charge Berry phase and the Abrikosov-Hofstadter 
problem but by the Wigner-style crystallization. This 
is precisely the opposite limit of the Maxwell equation 
((40(1 . in which it is the dual magnetic energy that deter- 
mines the pattern rather than simply reflecting the one 
set by the Abrikosov-Hofstadter condensation energy, en- 
coded in and communicated by j*. This is clearly 
seen in the Wigner crystal limit: consider an array of 
real-space pairs fixed in their positions. Such particles 
carry a unit dual flux and are completely invisible to 
vortices. As a result j* = and ((40() turns into the min- 
imization of the dual magnetic energy. In the real-space 
pair limit this dual magnetic energy is nothing but the 
original assumed interactions between the pair bosons: 
i Vir^^r,) = i / d\d^r'Bd{r)V{r-r')Bd{r'), since 
Bd{r) = J2i '^(i" ^ where {rj} are the pairs' positions 
- this is just the minimization of the potential energy in 
the Wigner problem. Consequently, the two descriptions, 
that of the Cooper pairs versus the real-space pairs, cor- 
respond to the two opposite limits of 1(401 139() , the first 
to the strongly type-II, the second to the strongly type-I 
limit of a dual superconductor. The density- wave pat- 
terns associated with these two limits are generally dif- 
ferent and can be distinguished in experiments. 

We now resume our discussion of the Abrikosov- 
Hofstadter problem. The simplest case, for all models, 
is / = 1/2 or a; = 0. For the HI model, the (anti)vortices 
prefer red plaquette locations and will be large com- 
pared to <!>{,, which can be safely ignored. The resulting 
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dual vortex array at half-filling is depicted in Fig. |S| 
The structure is a checkerboard with vortices in lo- 
cated on alternating black plaquettes. There is a single 
dual vortex per two unit cells of the Cu02 lattice, as 
expected for / — 1/2. Such close packing of dual vor- 
tices results in "empty" black plaquettes actually being 
occupied by dual antivortices (manifestation of the fact 
that our / = 1/2 Hofstadter problem does not break the 
dual version of time reversal invariance). There is, how- 
ever, no modulation in bond variables SA^j located on 
vertices of the blue lattice due to its peculiar relation to 
dual magnetic induction - all blue sites are in a per- 
fectly symmetric relation to the dual vortex-antivortex 
checkerboard pattern on black plaquettes, as is obvious 
from Fig. [3 This implies SAij — and the pseudogap A 
remains uniform despite $(r) 7^ 0. From this one would 
tend to conclude that there will be no CPCDW and no 
modulation in the local DOS. Still, there is a clear lat- 
tice translational symmetry breaking in the dual sector 
as evident from Figs. \E\ and Consequently, if we go 
beyond the leading derivatives, for example by including 
corrections to the Berry phase discussed around (|26|l , we 
expect that a weak checkerboard modulation will develop 
in quantities like the electron density Srii. The above is a 
special feature of the fluctuating LdSC which, however, is 
altered when we include the spin channel in our consider- 
ation. The Berry gauge field must then be restored - 
its coupling to nodal fermions induces antiferromagnetic 
order at half-filling and thus breaks the above symme- 
try in the leading order by a commensurate spin-density 
wave (SDW)i^. 

As the system is doped / decreases away form 1 /2 ac- 
cording to / = p/q — {1 — x)/2. The ground state energy 
of H38|l is particularly low for dopings such that g is a 
small integer, (integer)^ or a multiple of 2, reflecting the 
square symmetry of the Cu02 planes. Such dopings are 
thus identified as "major fractions" in the sense of the 
Abrikosov-Hofstadter problem. In the window of doping 
which is of interest in cuprates, these fractions are 7/16, 
4/9, 3/7, 6/13, 11/24, 15/32, 13/32, 29/64, 27/64, 
corresponding to dopings x — 0.125 (1/8), 0.111 (1/9), 
0.143 (1/7), 0.077 (1/13), 0.083 (1/12), 0.0625 (1/16), 
0.1875 (3/16), 0.09375 (3/32), 0.15625 (5/32), etc. Other 
potentially prominent fractions, like 1/4, 1/3, 2/5, or 
3/8, are associated with dopings outside the underdoped 
regime of strong vortex-antivortex fluctuations. In gen- 
eral, we expect that particularly low energy states cor- 
respond to fractions such that the pattern of dual vor- 
tices in $ can be easily accommodated by the underlying 
Cu02 lattice. Furthermore, we expect that the quartic 
repulsion in H38|l will favor the most uniform array of dual 
vortices that can be constructed from the ^-dimensional 
degenerate Hofstadter manifold. In the window of dop- 
ings one deals with in cuprates, these conditions single 
out doping X = 0.125 (/ = 7/16) as a particularly promi- 
nent fraction. At a; = 0.125 {q = 16), the dual vortex pat- 
tern and the accompanying modulation in can take 
advantage of a 4 x 4 elementary block which, when ori- 



ented along the x{y) direction, fits neatly into plaquettes 
of the dual lattice, as depicted in Fig. This 4 x 4 el- 
ementary block embedded into the original Cu02 lattice 
and containing 7 dual vortices (/ — 7/16) is clearly the 
most prominent geometrical structure among all the ones 
we have found in our study, both in its intrinsic simplic- 
ity and its favorable commensuration with the underlying 
atomic lattice. It is bound to be among the highly ener- 
getically preferred states in underdoped cuprates, as it is 
indeed found in the next section. 

Before we turn to the details of this energetics, we in- 
vestigate the signature in the electronic structure of the 
above 4x4 elementary block. Such signature could be 
detected in the STM experiments of the type performed 
in Rcfs. I^IMTol To this end, we compute the dual induc- 
tion Brf associated with the pattern of 7 dual vortices in 
Fig. [3 and from it the spatial modulation of the d-wave 
pseudogap SAij on each bond within the 4x4 supercell. 
With A thus fixed, we evaluated the local density of BdG 
fermion states of our LdSC model. The results are shown 
in Figs. I7I9I Note that we do our computations in the 
"supersolid" region of the phase diagram in Fig. 0] This 
enables us to use the BdG formalism with only moderate 
smearing in the coherence peaks coming from the gauge 
fields V and a and it also allows for rather direct com- 
parison with experiments. The downside is that we have 
to assume that the modulation profile in B^ remains the 
same as determined from our dual mean- field arguments. 
Considering the high symmetry of elementary block in 
Fig. lathis appears to be a rather minor assumption. 



VI. DUAL ABRIKOSOV-HOFSTADTER 
PROBLEM 

In this section we present the results of numerical anal- 
ysis of the HI model. Within this model the values of the 
matter field $(r) on the black sites are suppressed by very 
large m^, and ^(r) effectively lives on the lattice dual to 
the original copper lattice, that is on the red sites in Fig. 
12 Therefore, the dual fluxes reside inside the red plaque- 
ttes, shown in Fig. 1111 In this section we will drop the 
subscript d in order to make notation more compact and 
use A(r) for the vector potential corresponding to a uni- 
form dual magnetic flux equal to / = p/q. Modulation 
of the fleld, which will be determined numerically, is de- 
scribed by (5A(r). Within mean-field approximation, our 
problem then is reduced to minimization of the following 
Ginzburg-Landau lattice functional: 

Hd = -i^e'-^-'-+^+**-^-'-+^$*(r)$(r + ^) 

r.S 

-f^(m2|ci>(r)p + ||$(r)|4)+|^5i?^(r) , (41) 

r r 
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FIG. 5: Left panel: the circles depict the checkerboard array of vortices in <1> forming the ground state of Hamiltonian HI at 
half-filling (/ = 1/2 per plaquette of the atomic lattice). The empty black plaquettes are actually occupied by dual antivortices, 
brought into existence by a simple geometric constraint on the phase of $ in such a checkerboard array. The full green circles 
denote "dual vortex holes"—, i.e. the dual vortices that are missing relative to the half-filling checkerboard pattern once the 
system is doped to a; = 1/8 (/ = 7/16). The centers of these green circles form a square which defines the 4x4 elementary block 
of 7 dual vortices (black circles) per 16 sites of the Cu02 lattice discussed in text. Right panel: the most general distribution 
of modulations iJAij consistent with the vortex pattern just described. The green "crosses" correspond to "dual vortex holes". 
In general, there are six** phenomenological parameters JAij shown in different colors. The number of inequivalent sites of the 
black lattice is also six (see Fig. |7|l. 



where rr? = m^, 8 = {±i, ±y}, the hnk variables A are 
defined as 

A,r+5 = y dr-A . (42) 

Note that from this point on, we absorb the factors of 27r 
into the definition of A for numerical convenience and 
to conform with the Abrikosov dimensionlcss notation^i. 
The modulated part of the flux SBd{r) is given by the cir- 
culation of the corresponding SA around each plaquette 
as 

^Ar.r+x ^ ^Ar-\-x.r-\-x-\-y ^ ^Ar^x+y,r+y ^ ^-^r+y.r ■ (43) 

The minimization of Tid with respect to the link variables 
SA is equivalent to the solution of the following set of 
equations: 

= 2t|$(r)| |$(r -I- £)| sin (A,r+S + 5Ar,r+x + ar+x - "r) 

+ nliSBdir) - 5Ba{v - y)) 
= 2t|$(r)| |$(r + y)\ sin (A,r+y + <5A,r+y + cir+y - ay) 
+ nl{5Bd{Y ~ x) ~ SBdir)) , 

(44) 

where Ofr is the phase of the dual matter field $(r). Since 
the last terms in equations (|44|l can be identified as x and 
y components of the lattice curl, these are the lattice 
analogs of the (dual) Maxwell's equations in two dimen- 
sions, providing explicit lattice realization of Eq. (|40|l . 

Before we present the results of the numerical compu- 
tation, we will discuss briefiy the structure of the so- 
lutions that should be expected on general symmetry 



grounds. In the limit of infinite the gauge field 5 A 
does not fluctuate, and only <l>(r) should be varied in 
minimizing Tid- To understand why the solution for a 
general filling is inhomogeneous consider first a case of 
g — 0. Then the functional Tid simply describes a particle 
on a tight-binding lattice moving in a uniform magnetic 
field /. The corresponding Hamiltonian is 

HuoMr)^-t ^ e*-^'-'-+^$(r + ^) + m2$(r) . 

The minimization of the functional Tid is closely re- 
lated to finding the ground states of Hamiltonian -ffnof ■ 
Note that although the dual magnetic field felt by the 
particles is perfectly uniform, the gauge-field A(r) is not, 
regardless of the gauge used. Indeed, if this were the 
case the circulation of the vector potential A around the 
primitive plaquette would be zero by periodicity, which 
is not possible as the circulation is equal to the flux of 
the magnetic fleld 2'iTp/q through the plaquette. Thus, 
the Hamiltonian fffjof does not commute with the usual 
lattice translation operators. Instead, as noted by Peierls 
long time ago, magnetic translation operators Tr, gener- 
ating lattice translations complemented by simultaneous 
gauge transformations, must be constructed in order to 
commute with -ffHof ■ 

Unlike the ordinary translations, operators Tr do not 
commute. Rather, operators Tr. form a ray represen- 
tation of the translation group. The theory for irre- 
ducible ray representations of the translation group was 
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FIG. 6: The same as Fig. |K|except now the role of red and black plaquettes is reversed. This corresponds to either < or 
the (anti)vortex core location at the black plaquette being a local minimum of the vortex lattice potential V(r). In the latter 
case, the pictured array would be a metastable configuration of dual vortices at x = and x = 1/8, ultimately unstable to the 
true ground state at those dopings depicted in Fig. |S| The most general pattern of SAij is shown on the right. In this case, the 
"dual vortex holes" (solid green circles) correspond to the centers of SAij bonds shown in green. Note that in general there are 
si:!j— distinct parameters controlling modulations of the pairing pseudogap Aij. The number of non-equivalent (black) sites, 
however, is three (see Fig. |HJ. 



constructed by BrowrJ^. Alternatively, one can use the 
magnetic translation group introduced by Zak^, and use 
the ordinary representations of that group to classify the 
eigenstates of the Hamiltonian. 

On the lattice, the classification of the states is simple 
and for completeness, we demonstrate how the magnetic 
eigenstates can be constructed. Consider for example the 
Landau gauge = 0, Ay = 2TTxp/q, in which the unit 
cell spans q elementary plaquettes in ^-direction. Then 
the Peierls factors Ar.r+s , shown in Fig. EH are periodic 
modulo 27T with enlarged qxl unit cell. The Hamiltonian 
now can be written as 

HhoMx, y) = ?7i2$(j;, y) - t(^<P{x + l,y) + <i>{x - 1, y) 
+ $(a;, y + l)e'"f " + $(x, y - l)e-2"f . 

The Hamiltonian, obviously, remains invariant under 
transformations r — > r + y and r — > r + qx. Consequently, 
it can be diagonalized with the usual Bloch conditions 

$(r + g£) = e*«'=="$(r) (45) 
$(r + y) = e*'=«$(r) , (46) 

where (fc^;, ky) is the crystal momentum defined in a Bril- 
louin zone ^ x 2tt. Using these conditions, we rewrite 
the equation for the eigenstates as 

g{x) — t(^g{x + 1) + g{x — l) + 2g{x) cos{ky + 2nxp/ q)^ 

^Eg{x) , (47) 



where g{x) = $(a:,0). Now we have one-dimensional 
equation for g{x), where a; = 0, 1, 2, . . ., g — 1, that has to 
be solved with Bloch condition 

g{x + q) = e"'''-g{x) . 

Thus the problem of diagonalization is reduced to the di- 
agonalization of q x q matrix for each k. Note, however, 
that we did not exhaust all the information contained in 
the magnetic translation operators, apart from transla- 
tions by q lattice spacings in x direction. Consider an 
eigenstate described by crystal momentum k within the 
Brillouin zone and characterized with wavefunction g{x). 
Then function gi{x) = g{x + 1) is also an eigenstate of 
the Hamiltonian Huof with the same energy but with mo- 
mentum {kx,ky + 2'Kp/q). By repeating this operation, 
one finds that q states with crystal momenta described 
by the same but different ky: 

ky, ky + 2'rTp/q, ky + Anp/q, . . ., ky + (q - l)2TTp/q 

all have the same energy. Since p and q arc mutually 
prime, this set coincides with 

ky, ky + 2Tx/q, ky + 47r/g, . . ., ky + {q - 1)2tt / q . 

Occasionally, in the theory of magnetic translation 
groups, this is expressed by using a reduced Brillouin 
zone of size (2-k)'^ /q^, then every band is said to be q- 
fold degenerate. As a typical example, a dispersion for 
p/q — 1/4 is shown in the right panel of Fig. 1101 

In the gauge we just described, there are q degen- 
erate minima described by wavefunctions '^j{x,y) that 
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FIG. 7: The local density of states (LDOS) of the lattice d-wave superconductor with the modulated bond gap function Aij 
corresponding to the 4x4 supercell structure of Fig. I^at doping a; = 1/8. The calculations are done within the LdSC model 
of phase fluctuations. Note that the fluctuations of the gauge fields v and a lead to small broadening of the peaks but produce 
no significant changes, as long as one is in the superconducting state. The parameters are t* — 1.0, A = 0.1 and the variation 
in A from the weakest to the strongest bond is ±25%. The portion enclosed within a rectangle is enlarged in the central panel. 
The numbers correspond to the locations of Cu atoms within a 4 x 4 unit supercell on the Cu02 lattice depicted on the right 
panel. Note that Cu atom labeled 1 coincides with the location of the "dual vortex hole" in Fig. |^(full green circle) where A is 
the weakest, while the one labeled 6 corresponds to the Cu atom in the center of Fig. |K|which is surrounded by four strongest 
A's. The radii of the red circles indicate the magnitude of LDOS at E = — 0.27t* ( dashed line in the central panel). This 
pattern of LDOS is very robust in our calculations and is precisely the tight-binding analogue of the checkerboard structure 
observed by Hanaguri et al- (note that on a tight-binding lattice, the Fourier transforms at wavevectors 27r/4ao and 3 x 27r/4ao 
are not independent since 3 x 27r/4ao and — 27r/4ao are equivalent). The symmetry of modulations in Aij corresponding to Fig. 
imply that there are only six non-equivalent sites within the 4x4 unit cell. Note that modulation pattern at energies above 
A is actually reversed compared to the pattern at energies below A. Finally, the nodes remain effectively intact, in accordance 
with Ref. Ha . Finite LDOS at zero energy is entirely due to artificial broadening used to emulate finite experimental resolution. 
In the absence of such broadening, the LDOS remains zero within numerical accuracy (see Fig. |3J|. 



are located at = and ky = 2Trmp/q, where j = 
0,1,2,...((7— 1). Therefore, sufficiently close to the 
transition, the minimum of the functional Hd should be 
sought as a linear combination of the q degenerate states 
$j(r) and our problem is equivalent to minimizing the 
Abrikosov participation ratio 

^jj^^rl^WI where $(r) = Vc,$,(r) . 

(48) 

Since kx = for all (r), any linear combination of func- 
tions ^j{r) is periodic in x-direction: $(r -|- qx) — $(r). 
In addition, for each of the q ground states, ky is a mul- 
tiple of 27r/g, and consequently <I>(r -I- qy) — <I>(r). From 
these two properties we find that any linear combination 
J2j Cj $j (r) must be periodic in the qx. q unit cell. 

The minimization problem can be formulated equiva- 
lently in terms of coefficients Cj : 

H, = i?oEiQi'+ErS,..,3,..c;;^;^.3Q. , (49) 

where Eq is the ground state energy of -ffuof following 



from gTj) and 

^n,n..,u - E '^*n (r)$.3 (r) . (50) 

r 

Note that (|49|l itself has Ginzburg-Landau form with q 
order parameters Cj. In our case, the form of the ma- 
trix r^'') is dictated by our "microscopic" Hamiltonian 
Tid and corresponds to the Abrikosov participation ra- 
tio (|48|) . One could in principle generalize the theory by 
considering completely general form of P'^' compatible 
with the overall symmetry requirements. Such procedure 
is equivalent to introducing long-ranged quartic interac- 
tions with general kernel K{r,r' ,r" ,r'"): 

J2 Kir, r', r", r"')$*(r)$* (rO^-lr'^^Cr'") . 

Equation (|49|l is the most direct and convenient repre- 
sentation for numerical minimization of Tid in the limit 
of infinite Kd- However, in order to allow for intermedi- 
ate values of Kd and to be able to analyze the impact of 
various short-ranged terms describing interaction of dual 
fluxes 6Bd{r), we also have opted to follow a slightly dif- 
ferent route: The numerical results that are presented 
below are produced by direct numerical minimization of 
functional (|41|) with respect to both the dual matter field 
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FIG. 8: Left panel: The local density of states (LDOS) of the lattice d-wave superconductor with the modulated bond gap 
function Aij corresponding to the 4x4 supercell structure of Fig. El at doping x — 1/8. The parameters used are the same as 
in Fig. [7| Right panel: the local density of states at energy E — —0.27t (dashed line in the central panel). The radii of the 
circles are proportional to the LDOS at a given Cu atom. Note that inside the unit cell 4x4 there are three (rather than six, 
as in Fig. |7J non-equivalent sites. Even in absence of any underlying theory, this pattern of modulation appears to correspond 
to CPCDW simply by visual inspection. 
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FIG. 9: The local density of states (LDOS) of the lattice d-wave superconductor with the modulated bond gap function 
Aij corresponding to the 4x4 supercell structure of Fig. |S]at a; = 1/8, computed under the assumption of perfect particle- 
hole symmetry for low-energy fermionic excitations. The parameters used are the same as in Fig. |7| with the exception of the 
broadening, which has been suppressed here in order to demonstrate that the nodes remain essentially unaffected by CPCDW—. 



$ and the link variables 6Aa by imposing periodic bound- 
ary conditions on x Ny blocks with varying and 
Ny. At least in the vicinity of the transition, the largest 
unit cell one has to consider is qxq. It should be stressed 
that we found identical results using both approaches 
whenever direct comparison is possible (sufRciently large 
Kd and no additional short-ranged interactions between 
the fluxes). 

We perform numerical minimization of functional 141|) 
with respect to both the matter field ^(r) and the link 
variables SA^, by imposing periodic boundary conditions 
on n X TO blocks with variable to and n. The number of 
independent variables grows as Smn and the largest unit 
cells we were able to consider are 8x8. We used the 
conjugate gradient minimization technique with as many 



as 10^ — 10^ randomly chosen different starting points. 

We remind the reader that since the dual matter field 
variables <&(r) in the HI model live on the on red sites, 
the dual fluxes SBd{r) reside inside the red plaquettes, 
shown in Fig. ^] Note that each link of the original cop- 
per lattice is shared by two red plaquettes and therefore 
the enhancement of the c?-wave gap function A should 
be interpreted as the average between the fluxes at the 
neighboring dual plaquettes. Thus, at half- filling f = \ 
there is no modulation in A and no charge density mod- 
ulation emerging from our model HI even though there 
is a checkerboard pattern in the dual flux. The checker- 
board pattern remains the same for the entire range of 
parameters we were able to check. However, as explained 
in the previous section, since there is a definite symmetry 
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FIG. 10: Left panel: the Peierls link variables Aij in Landau gauge for q — 4. x denotes 2TTp/q. Both Ar,T+x and AT,r+y are 
periodic in ^-direction. Although AT,r+y increases monotonically with x, exp(i^r,r+4) is periodic with the unit cell shown by 
the dashed rectangle. Right panel: the dispersion E(kx, ky) of the lowest Hofstadter band for p/q = 1/4 in units of t* . There 
are g = 4 ground states at = and ky = 2n j/q, where j = 0,1, . . ., q — 1. For q = 4 the energy of the ground states is 
Eo = -2V2t. 



breaking in the dual sector by the checkerboard array of 
SBd{r) higher order derivatives and other terms not in- 
cluded in our dual Lagrangian are expected to generate 
a weak checkerboard modulation to accompany SBii{r). 

At field characterized by / = p/q = 7/16, which within 
HI model corresponds to doping x — 1/8, the structure 
of the configuration is considerably more complex. When 
restricted to a 4 x 4 lattice, the resulting pattern is the 
square lattice of "crosses" separated by four unit cells, 
shown in Fig. 1121 This configuration, however, is not 
the true global minimum. If a larger unit cell for mod- 
ulations of the dual fluxes is allowed, the energy can be 
additionally lowered by « 0.5% by distorting the ideal 
square pattern. The lowest energy we found corresponds 
to the quasi-triangular lattice of crosses shown in the 
right panel of Fig. ^1 The lowest energy state that we 
find has the symmetry of this quasi-triangular pattern 
for aU from 1.0 to 10^. 

The smallness of the energy differences, involving only 
few percents of the overall Abrikosov-Hofstadter energy 
scale, indicates that the state that emerges victorious can 
be changed by the additional short range interactions 
and derivative terms which we have routinely neglected. 
Not surprisingly, therefore, the precise energetics of var- 
ious low energy Abrikosov-Hofstadter states is decided 
by details. We find that inclusion of terms [^(r)!^ and 
dual density-density interactions |$(r)p|$(r -t- 5)p with 
moderate coefficients does not change the patterns we 
described. On the other hand, the inclusion of the terms 
describing short-ranged interactions between the dual 
fluxes produces signiflcant effects. An example of the 
typical pattern obtained by replacing the self-interaction 
term 

T^l^diSB,)' (51) 

r 

by 

r nn nnn 

(52) 



is shown in Fig. 1131 The parameters Fq and Fi in 
(|52|l are chosen to make the distribution of the dual flux 
somewhat smoother than what is demanded by 151() only. 
This suffices to bring the energy of the square pattern in 
Fig. ^|from just above to just below that of the quasi- 
triangular pattern of Fig. 1121 It is tempting to speculate 
that this slight additional "smoothening" of the dual fiux 
represents the combined effects of nodal fermions and 
Coulomb interactions present in real systems. Note that 
the symmetry and the qualitative features of this pattern 
coincide with our 4x4 "elementary" block conjectured to 
be the likely CPCDW ground state in the previous sec- 
tion (see Figs. ElandEJ. Once 5Bd{i) are translated into 
(5Ay the resulting pattern closely resembles the checker- 
board distribution of local density of states of "electron 
crystal" observed in the STM experiments, as illustrated 
in Fig. □ 

Let us now consider this observed STM checkerboard 
pattern in more detail, in light of our theory. Obvi- 
ously, our analysis being restricted restricted to the tight- 
binding lattice, we cannot describe the local density of 
states (LDOS) at positions between the sites of Cu lat- 
tice - the peaks in our theoretical LDOS are always lo- 
cated on top of Cu atoms. In contrast, the STM mea- 
surements by Hanaguri et al^ , have much better spatial 
resolution and can image the actual continuous atomic 
orbitals. Our tight-binding lattice results for LDOS(r, 
E) could be viewed as "coarse-grained" representation of 
the LDOS g{T, E) observed experimentally: 

LDOS(R,£)cx j g(Y,E)dY , 

where the integral extends over a square of size uq x ag 
centered around Cu lattice site R. Equivalently, the true 
continuum LDOS signal could be obtained by broadening 
our lattice LDOS around each Cu lattice site. 

A prominent feature of the checkerboard pattern in 
STM measurements that has received much attention 
is the presence of a pronounced Fourier signal not only 
at wavevectors and which correspond 

to the 4ao x 4ao periodicity just described, but also at 
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FIG. 11: Left panel: Any link of the black lattice (thick solid lines), which corresponds to A^j, is shared by two dual plaquettes 
shown in red. Within HI model, therefore, the enhancement or suppression of Aij is determined by the average of the fluxes 
trough the neighboring plaquettes of the dual lattice (red dashed lines). Central panel: distribution of 5Bd at half fiUing. In 
units of t, the parameters used are — 30.0, = —1.0, and g — 2.0. The energy per site is —3.695. Right panel: distribution 
of SBd at p/q = 1/4. In units of t, the parameters used are = 30.0, = —1.0, and g — 2.0. The energy per site is —3.564. 
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FIG. 12: Left panel: The pattern of the dual fluxes SBd for p/q = 7/16 with periodic boundary conditions in 4 x 4 unit cell. 
The unit cell 4 x 4 is repeated four times in horizontal and vertical directions for presentation purposes. The energy per site is 
—3.193 in the units of t and the parameters of the model are the same as in the previous figure. The positive(negative) values 
of SJid are shown in red (black). Right panel: the same but with 8x8 unit cell. The square lattice is distorted towards the 
triangular lattice. The energy per site is — 3.208t and is lower by 0.5% compared to the square arrangement. 



f^(|,0) and ^(0, |). While it is tempting to asso- 
ciate the 3/4 peaks with an entirely independent type 
of order, we note that for a periodic lattice of identical 
4ao X 4ao tiles, the Fourier transform is a discrete series 
with wavevectors Qn^ — \ (''^x i^y), where rix and riy 
are integers, irrespective of how complex is the internal 
structure of each tile. For a general structure of the tile 
all of the harmonics (n^;, riy) are present, and there is no 
a priori reason for the Fourier coefficients with = 3, 
ny = to be particularly small. The presence of a large 
signal at o is only natural, and no more unexpected 
than the weakness of Fourier harmonics at Q2,o and (5o,2- 
Next, in order to describe and compare the Fourier 
transforms of spatially broadened LDOS patterns in Figs. 
(|7I9|I and tunneling LDOS observed in experiments^, we 
introduce a simple model which approximates each bright 
spot inside the primitive 4ao x 4ao tile as 



r]{r) = exp 



2ttx 2-Ky . 

J (cos h cos 2; 

4ao 4ao 



While the specific functional form of the peak is not im- 



portant, our choice is convenient since g{T) has a period 
of 4ao X 4ao and a gaussian shape, centered at positions 
r = (4A'j;ao, 4iYyao), where and Ny are integers. The 
width of the gaussians is ~ ao/ 

The real-space tunneling LDOS pattern of Hanaguri et 
al. can now be represented by a function 

g(r) = yl,7(r)+i?^77(r + ^) + C^77(r + 5') , (53) 

where the first term represents the brightest peak of each 
tile, the term proportional to B represents the set of four 
second brightest peaks located at (5 = ±(1 -I- e)ax, ±(1 -|- 
e)ay and the last term represents the weak peaks at 8' = 
±ao(l -I- e)x ± ao(l + e)y (see Fig. [T^ . Parameter e 
equals zero if the real-space LDOS peaks are centered at 
the Cu lattice sites, while e = 1/3 if the maxima of all 
peaks, except for the central one, are displaced from their 
commensurate positions on top of Cu atoms - this is how 
the experimental data were interpreted in Ref. 0. We will 
show that for a wide range of parameters J, A, B,C, and 
arbitrary < e < 1/3 the Q3.0 peaks are much stronger 
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FIG. 13: Left panel: The pattern of the dual fluxes SB^ for the Hofstadter-Abrikosov problem at / = 7/16 (x = 1/8) obtained 
for the 8x8 unit cell. The parameters of the model are = 30.0, To = 6.0, Fi = —12.0, = —1.0, and g = 2.0. The 
4x4 checkerboard symmetry of the pattern is precisely what is shown in Fig. |7| Right panel: The distribution of modulations 
in the pairing pseudogap SAij corresponding to the pattern of fluxes shown on the left panel. This "tartan" pattern has the 
symmetry of Fig. I^with "dual vortex holes" coinciding with black "crosses". 



m 


e = 


e = 1/3 


m = mod 4 
m = 1 mod 4 
m = 2 mod 4 
m = 3 mod 4 


^ + 4B + 4C 
A + 2B 
A -AC 
A + 2B 


A + 4B +4C 
A + B-2C 
A + B-2C 
A + 4B+4C 



TABLE 1: Non-monotonic dependence gm,o of the Fourier 
transformed LDOS g(r) defined in 15311 . 



than peaks at Q2,o, although to explain other features 
of Fig. 2 in Ref. 9^ the "commensurate" choice e = 
appears to be more natural. For the momentum-space 
direction rt = shown as black circles in Fig. 2 of Ref. 
0, the Fourier coefficients are 



1 



9m,0 



4ao /•iao 



dxdy G{x, y)e 



The integral is elementary and the result for gm.o is 

27r(l + e)m' 



A + 2B + {2B + 4C) cos ■ 



2J 



(54) 

where Im is the regular modified Bessel function. Im{J), 
shown for several values of J in Fig. ^1 is a monotoni- 
cally decreasing function of m. The non-monotonic part 
of (7m. Oi denoted by gm,Oj is contained in the first factor of 
(|54|l . When e = (commensurate position of the peaks), 
this factor is a periodic function of m with period four, 
while for e — 1/3 the said period is equal to three. Table 
^summarizes the factors gm,Q for the two cases. 

We start our analysis with the incommensurate case 
e — 1/3. The third column of Table |l] implies that in 
this case ^i^o and 52,0 are equal and smaller than the 
component 33,0, as shown in the third column of the ta- 
ble. Experimentally, however, the Fourier component at 



Qi^o (the 1/4 peak) is roughly of the same magnitude as 
Qsfi, and it is the component at Q2,o that is the weak- 
est. To account for particularly small 172,0 ffi.o one has 
to select the value of J such that the monotonic function 
Im{J) decreases rapidly in the region between m = 1 and 
m = 2, Fig. 1141 indicates that this is the case for 3 < J. 
This choice, however, also dramatically suppresses the 
Fourier transforms at m = 3 and to = 4, both of which 
are rather large in experiments. 

While the qualitative features of the Fourier trans- 
formed experimental LDOS could possibly be reconciled 
with e = 1/3 by fine-tuning parameters A,B, and C, 
the situation might in fact be better described by assum- 
ing the commensurate location of the peaks, in registry 
with Cu atoms (e — 0); see the second column the Ta- 
ble ^ In this case, harmonic to = 2 is automatically 
suppressed compared to the 1/4 and 3/4 Fourier com- 
ponents. The suppression of Fourier component at 2/4 
can be qualitatively understood as follows: This Fourier 
component is determined by the overlap of the LDOS sig- 
nal and cos(2^a;). Obviously, the maxima of the LDOS 
signal correspond to alternating maxima and minima of 
cos(2^a;) and destructive interference of the two func- 
tions occurs. For harmonics 1/4 and 3/4 (and of course 
for 4/4) the overlaps are significant and their Fourier co- 
efficients are larger. Peaks corresponding to larger val- 
ues of m are strongly reduced due to the monotonic de- 
pendence Im{J) of the Fourier transform on to. This 
suppression can serve as an estimate of parameter J: 
visually, the last discernible peak in Fig. 2 of Ref. 
appears at Qs.o, which places J in the range of 5 — 10 
(see right panel of Fig. I14|l . In the right panel of Fig. 
El the spatial Fourier transforms of the LDOS (|53|) with 
parameters J = 7, A = 1, B = 0.5, and C — and com- 
mensurate placement of the peaks e = is shown. This 
example illustrates that the major features of the Fourier- 
transformed LDOS obtained by Hanaguri et al.^ are ro- 
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FIG. 14: Left panel: Function g{r) used to emulate the LDOS signal with parameters J — 7.0, a = 1.0, b — 0.5 and c — 0.2. 
Center panel: Regular modified Bessel function Im{J) shown as a function of its index m for several values of fixed J. Right 
panel: Fourier coefficients g{qx, 0) for LDOS pattern shown in the left panel. The reader should compare this with the similar 
plot in Fig. 2 of Ref. |9|. 



bust properties of an elementary tile of size 4ao x 4ao and 
nine peaks occupying commensurate locations at sites of 
the Cu lattice as depicted in our Figs. H7I9I) : the large 
magnitude of 3/4 peak is simply a higher harmonic de- 
scribing the characteristic intra-tilc structure. 



VII. CONCLUSIONS 

Our main goal in this paper is to devise a more re- 
alistic description of a strongly quantum and thermally 
fluctuating cZ-wave superconductor, based on the theory 
of Ref. I3. Such description applies not only to long 
distance and low energy properties, which are the pri- 
mary domain oS^, but also to intermediate lengthscales, 
of order of several lattice spacings, and to energies up 
to the pseudogap scale A. This enables us to use the 
t heory to address the experimental observations of Refs. 
IS Utlirl The charge modulation observed in those exper- 
iments is attributed to the formation of the Cooper pair 
CDW, the dynamical origin of which is in strong quan- 
tum fluctuations of vortex-antivortex pairs. These quan- 
tum superconducting phase fluctuations reflect enhanced 
Mott-Hubbard correlations in underdoped cuprates as 
doping approaches zero. Quantum fluctuating hc/2e 
(anti)vortices "see" physical electron as a source of a 
half-quantum dual magnetic flux and the theory of the 
CPCDW can be formulated as the Abrikosov-Hofstadter 
problem in a type-II dual superconductor^^. An XY- 
type model of such a dual superconductor appropriate 
for a lattice d-wave superconductor is constructed, both 
for thermal and quantum phase fluctuations. The specific 
translational symmetry breaking patterns that arise from 
the dual Abrikosov-Hofstadter problem are discussed for 
various dopings x, which determines the dual flux per 
unit cell of the Cu02 lattice via / = p/q = (1 — x)/2. 
In turn, the spatial modulation of the dual magnetic in- 



duction Brf corresponding to these Abrikosov-Hofstadter 
patterns is related to the modulation in the gap func- 
tion of the lattice d-wave superconductor and is used to 
compute LDOS observed in STM experiments. A good 
agreement is found for x = 1/8 (/ = 7/16), which is the 
dominant fraction of the Abrikosov-Hofstadter problem 
in the window of dopings where our theory applies. 
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APPENDIX: TWO ALTERNATIVE 
DERIVATIONS OF DUAL 
ABRIKOSOV-HOFSTADTER HAMILTONIAN 

In this Appendix we present two different and self- 
contained derivations of the dual Abrikosov-Hofstadter 
Hamiltonian (|37I38|I . either one of which can serve as 
an alternative to the derivation given in the main text. 
The first approach is somewhat more detailed and in a 
sense more "microscopic" since it uses a quantum vortex- 
antivortex Hamiltonian as a springboard to derive the 
effective dual field theory (|36|l . In turn, such vortex- 
antivortex Hamiltonian in principle can be derived from 
the (still unknown) fully microscopic theory of cuprates. 
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Incidentally, this derivation is the (2+1) dimensional ana- 
logue of the 3D case presented in the Appendix of Ref. 
The second derivation follows the familiar Villain approx- 
imation to the XY model and applies it to our specific 
situation. The Villain approximation is less "realistic" 
but provides a transparent and systematic way of deriv- 
ing dual representations of XY-like models. 

In both derivations the starting point is the effective 
(2-|-l)D XY model of a quantum fluctuating c?-wave su- 
perconductor: 

i i nn 

- Jl X COs((pj - ipj)~ J2 X COs((^i - (pj) 
rnnn bnnn 

+ Lnodal[cOs{ipi ~ ipj)] + Lcoic , (A.l) 

where tpi (r) is the fluctuating phase on a site of the blue 
lattice in Fig. |21 the first (imaginary) term is the charge 
Berry phase corresponding to the overall flux / through a 
plaquette the blue lattice, J is the nearest neighbor XY 
coupling, t/i(2) are the next nearest neighbor XY cou- 
plings along red (black) diagonals, inodai is the contri- 
bution of nodal fermions, and Lcore denotes core contri- 
butions arising from small regions around vortices where 
the pairing pseudogap is significantly suppressed - this, 
for example, generally includes the mass term l|25|l . the 
energy cost of core-core overlap, the Bardeen-Stephen 
core dissipation, etc. The reader should bear in mind 
that the last effect is small in cuprates, as explained in 
the main text, and will be neglected in the Appendix. 
Furthermore, we are neglecting vortex interactions with 
the spin of low energy nodal fermions, represented by 
the Berry gauge field a; this is justified away from the 
critical point. The results below are easily adapted to 
the extended s-wave pairing symmetry. Similarly, both 
derivations are straightforwardly applied to a yet simpler 
case, a fluctuating s-wave superconductor: 



LxY^^fJ:^.+^Y.^^-JY: 



cos( 



(A.2) 

which was used in the main text as a pedagogical exam- 
ple. 

Just as the starting points of two derivations coincide, 
their final product, the effective dual theory at long and 
intermediate lengthscales, will also turn out to be the 
same. 



1. Tid from a "microscopic" vortex-antivortex 
Hamiltonian 

The quantum partition function of a phase fluctuating 
superconductor is: 

Ziy ^ j V^,ei,v{- j^rLj^ylip^iT)]) , (A.3) 



where the functional integral / Vipi{T) runs over phase 
variables fiir) such that exp(i(^i(T)) is periodic in the 
interval r € [0, (3]. 

The difficulty in computing l|A.3p is twofold: the fact 
that (fiir) is a compact phase variable, defined on an 
interval [0, 27r), rather than an ordinary real field tak- 
ing values in [— oo,-|-oo], and the cosine functions in 
Lj^Y that couple phases on different sites in a non- 
linear fashion. To deal with the problem one approx- 
imates the cosines with quadratic forms. One popu- 
lar approximation on a lattice is due to Villain and 
will be discussed in the next subsection. In continuum, 
the approximation amounts to replacing cos{ipi — fj) — *■ 
1 — (a^/2)(V(/j)^ + • • • , where a is the lattice spacing 
and ip{x) is now a function in continuous (2-)-l) dimen- 
sional spacetime. Its compact character is enforced by 
writing d^ip{x) d^x{x) + {^^l^{x))v, where x is an or- 
dinary real field and (9^(^(a;))^ is the part of the phase 
associated with vortices, defined via V x (V(/3(r, t))i, = 
2^ <5(r (r) ) - 27r E„ '5(r - (r) ) , with {r^i^'^) (r)} 
being (anti)vortex positions. Simultaneously with this 
decomposition of d^ip{x), J "Dipi is replaced by the func- 
tional integrations over x(^) and (anti)vortex positions 
{rQ^"-*(T)}. We are assuming here that the (anti)vortices 
of topological charge ±1 dominate the fluctuation behav- 
ior in the regime where the amplitude of the pseudogap 
A is large and stiff, allowing us to safely neglect topo- 
logical defects corresponding to vorticity ±2, ±3, . . . due 
to their higher core energies. This assumption simplifies 
the algebra considerably. Furthermore,^ this assumption 
is natural within the theory of Ref. ;§: The prolifera- 
tion of defects of high topological charge is equivalent to 
strong amplitude fluctuations and the eventual collapse 
of the pseudogap - as long as we are in the pseudogap 
regime the ±1 (anti) vortices are the only relevant excita- 
tions. With these changes in place, the partition function 
"3)1 finally takes the form: 



00 00 



^XY 



Af„=0 Wa=0 



E E whji / n 



a,7— 1 



J I?r;(r)exp(- J d'xL'j,y[x, {C^'Hr)}]) , 



{r»(0)} = {r»(/3)} 



(A.4) 



where / d'^x = Jq dr J d'^r and the set {rQ^"^(r)} con- 
taining A^^(a) (anti)vortices at r = coincides with the 
one at T = /3, to ensure proper periodicity of exp(i(^(r, r)) 
in imaginary time. Lastly, 

+ ^(Vx + V^.)2 + L„odai + LLc[K^"Hr)}] , (A.5) 

where /(r) = fi5[r — Hi), {Ri} are the sites of the 
blue lattice, J = J + Ji + J2, Kq has been rescaled by 
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a^, and 



„ |r-rX(T)P 
^ (r-r°(T)) X z 



7 



, (r-rQxz 
|r-r')P 

(r-rj;(T)) X z 



J(r,r), 



E 



(r - r^(T)) X z 
|r-r-(r)P 



= - / 



r-rX(r)|2 

, (r — r') X z 



|r-r')p 



(A.6) 



In (|A.6p we found it useful to introduce vorticity density 
and current: n(r, r) = (5(r-r^(T))-X;^ ^(r-r°(T)), 
J(i-, r) = ra-^lr-r- (r))-E^ r!;,5(r-r-(r)), in terms 
of which, when combined with vortex particle density 
and current p(r, r) = ^(r- ^aM) + E7 ^(i"- i"7W)' 
iir, r) = r-J{r - C(t)) + " r^W), we can 

write 



1 / dr'^ \ 2 

2E^H^) '5(r-r^;) + ^ctc , (A.7) 



where 



Hi,, = V(r)p{r,T) + ^p{r,r) J d\'V^'\r,r')p{r\T) 
+ ij,(r,r) I dV'F,?'(r,r')jKr',r) 
+ in(r,T) y d2^'V(2)(r,r')n(r',r) 
+ ij,(r,r)y'dVv(?)(r,r')J;(r',r) + (...) . (A.8) 

(Anti)vortex mass terms appearing in ljA.71) have been 
introduced already in Section IV (see the discussion sur- 
rounding Eq. H25|) '). while -ffcoro represents a systematic 
expansion in vortex core density, including single core, 
two core terms and so on. All the short range terms that 
arise from expanding the cosine functions ljA.l|) in con- 
tinuum limit have been absorbed into -ffcore> was our 
habit throughout the text - in particular, V{r) is just 
the vortex potential on the blue lattice, containing cru- 
cial information on d-wave pairing, which is extensively 
discussed in Section III. 

We can now integrate out x? the regular (XY "spin- 
wave") part of the phase. The quadratic phase stiffness 
terms in (|A.5p are decoupled as: 



iW ■ (Vx + V(/3^,) 



1 -.2 1 



via the Hubbard-Stratonovich vector field W = (Wq, W) 

(note that we have set the "dual speed of light" J/Kq 
to unity). Integration by parts gives iW ■ dx ^ —i{d- 
W)x and is followed by functional integration over x{^)i 
resulting in the local (5- function constraint 6{d ■ W). The 
constraint is solved by introducing a non-compact gauge 
field Ad such that W = d x Ad, ensuring d-W — d - {d x 
Ad) = 0. What remains of ljA.9|) is further transformed 
by another partial integration: 



^(r,^)- i{d X Ad) ■ {d^). + 



1 



2K, 



1 

2J' 



idxAd)l + i~{dxAd)\ 



lAd-dx {dip), + —{d X Ad)l + -~{dx Ad) 



2Ko 



2 J 



(A.IO) 



where (d x Ad)o.± denotes temporal and spatial compo- 
nents of 9 X Ad, respectively. Now observe that HA.6() 
implies d x {dtp), = (27rn, 27rJ). This allows us to fi- 
nally write the partition function of the quantum vortex- 
antivortex system as: 



00 00 



N„,N„ 



X / I?r°(r)exp(- f d^^xL^) , (A.ll) 

i{r»(0)} = {r»(/3)} J 



where LjJ equals 



dr 



2TTiAdon - i27r(A[,"' + A^) • J 



+ -^{dxAd)l + X{dxAd)l , (A.12) 
2Ko 2J 



and V X A 



(0) 



B 



(0) 



f{r)z. 



Eqs. ljA.lllA.12|l are an important result of this 
Appendix. We recognize Z!^ as equivalent to a parti- 
tion function of two species of non-relativistic quantum 
bosons expressed in the Feynman path integral represen- 
tation over particle worldline trajectories. These vortex 
and antivortex bosons have identical mass M and carry 
dual charges +2tt and — 27r, respectively, through which 
they couple to a dynamical gauge field Ad- The dual 
photons of Ad mediate long range "electrodynamic" in- 
teractions between the bosons, which are just the famil- 
iar Biot-Savart interactions between (anti)vortices. In 
addition, the particles interact through an assortment of 
short range interactions contained in iJcoro IIA.8|) . Fur- 
thermore, inodai describes the interactions generated by 
nodal Dirac-like fermions which will be included explic- 
itly once we arrive at the dual representation of l|A.12|l . as 
detailed in Section IV. Finally, the vortex potential V{r) 
contains important information about the underlying lat- 
tice structure and the symmetry of the order parameter, 
as emphasized throughout the text. 
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captures the general description of the 1 ^(o) , jg ^ \^ 12 

■antivortex svstem. aoDhcable to aU su- 9M'^ ^^'^ d- ^^'^ ^' ^' 



We have derived Z!^ as a continuum limit approxima- 
tion to the partition function Z'^y HA.3() of the quantum 
XY-type model (IA.1|I . Actually, in real cuprates and in 
all other physical systems, the opposite is true: It is the 
quantum XY-type representation that is an approxima- 
tion to Z^. 

quantum vortex-antivortex system, applicable to all su- 
perconductors and superfluids whose order parameter is 
a complex scalar. To appreciate this, imagine that for 
each given configuration of the phase, with (anti)vortex 
positions fixed in Euclidean spacetime, the microscopic 
action of a physical system is minimized with respect to 
the amplitude, after all other degrees of freedom have 
been integrated out. The subsequent summation over 
all distinct (anti)vortex positions leads to precisely Z^ as 
the final result (again, we remind the reader that the core 
dissipative terms will also generically appear in Z{f but, 
being small in underdoped cuprates, are neglected here as 
explained in the main text). In practice, this procedure 
is difficult to carry out explicitly and the actual values of 
various terms that enter Zf; are hard to determine from 
"first principles" . This is particularly true for core-core 
interaction terms appearing in i?coro- ^'^^ '^^^ purposes it 
will suffice to approximate T^(^^(r,r') g6{r — r'), where 
g > 0, and drop the rest. 



basis of coherent states: 



{dr + iedAdo)-^, 



There is one crucial feature which distinguishes Z.y 
from the standard Feynman partition function: The vor- 
tex and antivortex quantum bosons are not conserved. 
As particles make their way through imaginary time, 
vortices and antivortices can annihilate each other; sim- 
ilarly, they can also be created at an instant in time; 
this is depicted in Fig. |21 All such processes of cre- 
ation and annihilation proceed in pairs of vortices and 
antivortices. Consequently, while the individual number 
of vortices and antivortices is not conserved, the vortic- 
ity, measured by dual charge Cd = ±27r, is conserved and 
the gauge symmetry associated with Ad is always main- 
tained (unless, of course, it is spontaneously broken by 
a dual Higgs mechanism in dual superfluid). In other 
words, these nonrelativistic (anti)vortex bosons propa- 
gate through spacetime permeated by a vortex-antivortex 
condensate, of strength A^. 



2Af' 
1 

2M 

f(i*.„p+i*.n^+L„odai+2^ 

I 



|(V - iedAf^ - iedAd)-^a\^ - A^.d*. 



I*al 



l^idxAd)l)] , (A.13) 



where the meaning of various terms is straightforward in 
light of our earlier discussion; we have consolidated the 
notation so that ed = 27r, K^^ = {Kq, J ^ J), the chemi- 
cal potential for vortices is /z-u = IJ-a, and g describes the 
short range core-core repulsions. The "vortex-antivortex 
pairing" function Ay is crucial since it regulates the fre- 
quency of vortex-antivortex pair creation and annihila- 
tion processes. 

The form of Z^ can be further simplified by exploit- 
ing the vorticity conservation law. We observe that the 
action in ljA.13|) is invariant under gauge transforma- 
tions: exp«)«'i,, ^ exp(-iC)^'a, Ad 
Ad — {dO/cd- This prompts us to introduce bosonic 
"spinors" ^ — (^*,\E'q) which carry a conserved dual 
charge Cd and couple minimally to Ad'- 

Zv~^ j V-i^VAdVaexp[- j d^x{-^C^ 

+ inodai + ^a2 + ^(ax , (A.14) 



where 

C = 



{dr + icdAdo) Ay 

a; -idr+iedAdo)_ 

--i--(V + iedA^°^ -I- iedAdf - fiv 
ZM 



V{r)+ia)l , 
(A.15) 



and a Hubbard-Stratonovich scalar field a was deployed 
to decouple short range repulsion. Now we set Ad ^ 
and ignore inodai~ they will be easily restored later - 
and note that the integration over vortex matter fields ^' 
gives: 



Feynman path integrals are beautiful but difficult to 
calculate with. Following the standard mapping"*^ we 
can express Z^ as a functional integral over complex fields 
^'i,(r, t) and 5'a(r, t), which are the eigenvalues of vortex 
and antivortex annihilation operators, respectively, in the 



det 



Da exp ( 



d^x 



1 



Al 



25 
V{v) 



2M ^ 



+ V{v) 



The above partition function has a transition at /z„ = 
— |A^| (we are assuming that the minimum of V{v) oc- 
curs at zero, with E]i, or i?^ as the case may be, having 
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been absorbed into For /x„ < — |A„| the system of 
(anti)vortex bosons is in its "normal" state, with (^f) = 0. 
For fly > — |At,| (anti)vortex bosons condense and (4') 
becomes finite. This is nothing but the dual description 
of the superconducting transition discussed in the main 
text. In the general vicinity of the transition it is use- 
ful to introduce = /x^ — |A„p, where now m? > 
and < indicate dual normal and superfluid states, 
respectively. Focusing on distances longer than \/ A^M 
and energies lower than A„, the determinant in HA.16() 
can be further reduced to: 



det 



\A,\+V{r)+ia^ 
2M 



-™'-^[(W)+*(vV)] 

+ 2|A„|(y(r) + la) + {V{r) + laf] +(•••) , (A.17) 

with additional terms (• ■ ■ ) contributing unimportant 
derivatives. 

By setting ^(r) and cr(a;) to zero, we observe that 
the above expression assumes the form of the partition 
function determinant for a system of relativistic quantum 



bosons of mass 



det 



2y2 



with the 



speed of "light" c = a/|A^|/2M set to unity henceforth. 
The terms involving ^(r) and cr(x) describe the under- 
lying potential and various short range interactions of 
these relativistic bosons. We therefore can reexpress the 
determinant (jA.17p as a functional integral over the rela- 
tivistic boson field ^{x); this is a faithful representation 
of the original (anti)vortex partition function at distances 
longer than \/ AyM and energies lower than A^ : 

exp[- + i2TrAd)^\^ + m^{r)\^\^ 



+2|<i>R|A,| 



~V{r) 



2m]M + A/ 



^(9xA,)2],(A.18) 



where m'^{r) = + 2|A^,|V^(r) + ^(r)^ - {W'^V{r)/M) 
and we have restored £nodai and dual gauge field A^, 
through covariant derivatives d ^ D = {Dq, D) = {do + 

i2edAdo,'^ + i<^d-^''p + i&d-^d)- The minimal coupling 
of Ad is mandated by dual charge conservation: $ 
exp«)$, Ad^Ad- {dO/ed- 

The dual partition function Zd (|A.18|) is the final result 
of this subsection. It describes the system of relativis- 
tic quantum bosons of mass m and charge Cd — 27r in 
a magnetic field = V x A^^-*. The virtual particle- 
antiparticle creation and annihilation processes in the 
vacuum of this theory are nothing but quantum vortex- 
antivortex pair excitations evolving in imaginary time 
(see Fig. In the "normal" vacuum, > 0, the aver- 



age size of such pairs is ■ 



This is just the supercon- 



ducting ground state of physical underdoped cuprates. 
For < 0, this "normal" vacuum is unstable to a 
Higgs phase, with a finite dual condensate ($). The 
vortex-antivortex pairs unbind as infinite loops of virtual 
particle-antiparticle excitations of $ permeate this dual 



Higgs vacuum - this is the pseudogap state of cuprates. 
The integration over the Hubbard-Stratonovich field a 
produces a short range repulsion ^(4| A^, p)(/|<I>|* followed 
by an assortment of other short range interactions includ- 
ing V^(r), powers of |<I>| higher than quartic and various 
derivatives. All these additional interactions are irrel- 
evant in the sense of long distance behavior but might 
play some quantitative role at intermediate lengthscales. 
For simplicity, we shall mostly ignore them in this paper. 
Finally, with the change of notation 4|At,pg g and 
'Cnodai incorporated into the self-action for Ad as detailed 
in Section IV, the dual Lagrangian in IjA.lSI) reduces to 
Cd H36I) . The arguments in the text can then be followed 
to arrive at Hd I|37I38|I . 



2. Tid in the Villain approximation 

A useful approximation to the ordinary XY model is 
due to Villain. In this approximation the exponential 
of the cosine function is replaced by an infinite sum of 
the exponentials of parabolas. We will illustrate this ap- 
proximation first for the 2D case and follow up with the 
(2-f 1)D quantum XY-like model. 



a. Classical (thermal) phase fluctuations 

Here we apply the Villain approximation to our clas- 
sical XY model of a phase fiuctuating two-dimensional 
d-wave superconductor H14|l . As we will see presently, 
the final Coulomb gas representation coincides with 
- however, the advantage of the Villain approximation is 
that it will allow us to obtain explicit expressions for the 

r(b) 

core energies iJc in terms of the coupling constants of 
the original Hamiltonian (|14|1 . 

Denoting the sites of the blue lattice hy p — {x,y), the 
partition function of the model can be written as 

Z = J (Y\. '^'^pj X! ["^ (cos(Va;0p) -I- COs{Vy(l)p)) 

+ J12{P) (COS((/)p + £ + y - 0p) -I- COS((/)p+y - (l)p + x)) 

(A.19) 

where the lattice operator V is defined according to 
V5/(r) = fr+s — fr, and coefficients Juip) denote Ji 
(or J2) if p is in the lower left corner of a red (or black) 
plaquette. 

In Villain approximation the exponent of a cosine is 
replaced by a sum of Gaussian exponents that has the 
same periodicity 27r: 



+00 



exp[/3cos7] ~ i?v(/3) ^ exp[- 



(7 - 2?™)^ 



(A.20) 
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The fitting functions /3y(/3) and Rv{P) are determined 
by the requirement that the lowest Fourier coefficients of 
the two functions coincide. In particular, the function 
Pv iP) has the following asymptotic behavior for low and 
high temperatures: 



f3v{P) 



(3 ior(3 > 1 

21n|) for/3 <1 



(A.21) 



The sum over n in (|A.20p can be transformed by de- 
coupling the quadratic term in the exponent via the 
Hubbard-Stratonovich transformation: 

exp[/3cos7] cx V exp[- ———n'^ + ijn] . (A.22) 
^ ^Pv(P) 



n— — oo 



By introducing integer fields Ux{p), Uy[p) and w±{p) to 
apply (|A.22p to the cosine terms in (|A.19|1 we obtain: 

giEp["<:>(p)VQ0p + IU+(p)(0p+ii + j-0p)+tO_(p)(0p+j-0p + ^)] 

(A.23j 

where a — x,y, and coefficients Juip) and J' are defined 

as 



J' = PviJ) 

JMp) ^ PviJMp)) 



(A.24) 
(A.25) 



The integration over the angles 0p can be performed after 
applying the discrete analog of the integration by parts 
to the sums in the last exponent in (|A.23|I : 



E = E fip)i9{P + 9{P)) 

X X 

= J2(fiP -£)- f{p))g{p) = -Y.^f{p) 9{p) ■ 

X X 

(A.26) 

Note that we distinguish between the "right difference" 
operator V and "left difference" operator V. Integration 
over the phases (pp yields the following expression for the 
partition function: 



E 



■{p) + <AP) . "■1(P) + ™1(P) 



-H 777- 



^■^i2(p) \S{W-v{p) + ...) 



_ (A.27) 
where V-v(p) denotes two-dimensional lattice divergence 
of v(p) and dots denote 

(w+(p) - w+{p -x-y)) + {w^{p - x) -w^{p- y)). 

We rewrite the constraint appearing as the Kronecker 
delta function in the sum as 

V • (v w[u'+, = , 



where w = {wx,Wy) denote the following linear combi- 
nations of integer fields ■w± (p) : 

w+{p) + w+{p - y) w^{p) +w-{p~ y) 

Wx — 7^ 7: 

2 2 

(A.28) 

w+{p) + w+{p — x) W-[p) + W-{p — x) 



(A.29) 



The constraint can then be resolved as 



where 



V = b — w[w-|_, w_] , 



b = (VyA, ~V,A) 



(A.30) 



and A(p) has the meaning of the time- like component 
of a vector potential. At this point it is useful to pause 
and establish a simple geometrical interpretation of the 
various fields we have introduced. Variables Ux[p) and 
Uy (p) are coupled to the phase differences (fip+x — 4>p and 
4'p+y ~ 4'p ^-iid therefore they reside on the links ema- 
nating from p in positive x and y directions respectively. 
Integers A(p), on the other hand are related to link vari- 
able Ux{p) through the difference A(p) — A(p — y). Con- 
sequently we must associate A(p) with the centers of the 
blue plaquettes, which coincide with either red or black 
sites. Nevertheless, we will continue to use notation A(p) 
tacitly implying that p refers to the lower left corner of 
the (red or black) plaquette associated with A. Having 
resolved the constraint, we find that the partition func- 
tion Z can be written in terms of integer-valued fields 
w±[p) and A(p) only: 



E «^pE 



A, 104 



(b[A]-w[it;±])i 
2J' 



w\{p) 



wl{p) 



2J{2(P) 

(A.31) 

To obtain the description in terms of continuous rather 
than integer- valued fields A(p), we use the Poisson sum- 
mation formula: 



OC pQC ^ 

E /(A) = dAY: 

— oo I 



/(A) 



(A.32) 



A= — oo 



/ — — oo 



The partition function Z assumes the following form: 

n rfA(p) ^ exp ^ 2^i/(p)A(p)} exp{i^[b[A(p)]] , 

p i(p) p 

(A.33) 

where we have defined a functional exp{_F[b(p)]} accord- 
ing to 



^exp 



E 



(b-w[w±])2^ wl{p)+wl{p) 



2J' 



2J{2(P) 



(A.34) 



32 



In the limit when constants Ji and J2 are infinites- 
imally small, only the configurations w±{p) = con- 
tribute to the F[h]. This limit, which corresponds to the 
usual 2D XY model, is described by partition function 
Zq given by 



00 

-00 P 



( 2.a{p)K{p) - 



(A.35) 



For finite J12 we must resort to approximate evaluation 
of the functional F{h{p)): 



exp{i^[6,(p), = ^^P[- E 



w± p 

,,2 f „\ I „,.2 



2J' 



2J(2(P) 



(A.36) 



The quadratic terms containing b can be decoupled using 
the Hubbard-Stratonovich transformation: 



:,F[Hp)] - 



l[dZ^{p)dZy{p)X 



w± p 

iZa 



wlip) + w'i{p)- 



(A.37) 



where a — x,y. Using explicit expressions for Wa we 
obtain 



e^I"^! = / l[dZ,{p)dZy{p): 



w± p 



J'Zl{p) wl{p)+wl{p) 
2 2J'M 



-iZa{p)ba{p) 



-w+{p){Z^{p) + Z,{p + y) + Zyip) + Zy{p + x)) 



-w^{p){-Z, (p) ~Z,{p + y)+ Zy [p) + Zy{p + x)) 



(A.38) 



The sums over w± (p) can be performed by employing 
the Villain approximation (jA.20p backward. Note that 
the coupling constants Jj_ are restored to the original 



values of coupling constants J±: 

^fiHp)] = / Y[dZ,ip)dZyip)x 
p 

p 

^ , / „^ Z^jp) + Z.,{p + y) + Zy{p) + Zy{p + x) 

+ ■Jl2{P) [COS 



cos ■ 



-Z.,{p) -Z^{p + y) + Zy{p) + Zy{p + x] 



(A.39) 

To quadratic order, the expression in the exponent is 
p 



+ iZa{p)baip) - Jl2(p)x 



(Z,(p) + Z,{p + y)f + (Zyip) + Zy{p + x)f 



(A.40) 

Note that Z^ and Zy components are completely decou- 
pled at the quadratic level. To proceed, one can double 
the unit cell, in which case the expression in the exponent 
becomes diagonal in the momentum space. Alternatively, 
one can use an equivalent, but technically simpler proce- 
dure of keeping the original unit cell. In this latter case, 
the momentum space problem reduces to the diagonal- 
ization of a 2 x 2 matrix connecting modes at wavevectors 
q and q — g, where g — Tr{x + y). 

It is convenient to represent Ji2{p) as 



JM = 



J I + J 2 



J1-J2 



= J + 5Je'«P 



After Fourier transformation, the exponent in HA.40p be- 
comes 

■7 f \U f \ ry I .ry I yJ' + J {l + COS Qy) 

iZ^{q)bx{-q) - Z^(q)Zj,(-q) 2_ 

id J 

^sin9yZj;(g- q)Z^(q) + (a; ^ y) . (A.41) 

Since Za{p) and ba{p) are real, their Fourier components 
satisfy 



Z„(-q) = Z:(q) 

fea(-q) = ^a(q) • 



(A.42) 
(A.43) 



In the last expression for the partition function we found 
that the terms with Z^ and Zy decouple and we thus can 
integrate over Zx{p) and Zy{p) separately. 

The expressions in this and especially the next sub- 
section can be significantly economized by using a check 
mark to denote two-component vectors: 



&(q) 



>(q-g), 
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Using this notation, the contribution due to Zx{p) can 
be written as 



-6^(-q)Z,(q)--Zi(-q)(...)^.(q) 



(A.44) 



where {...) — J' + J + J cos qya^ + S J sin qya2 and super- 
script T denotes the transpose of a matrix. Now can 
be integrated out. Apart from the overaU normalization 
constant, F[b] is given by: 



"^"^ ^6l(-q)G(g,)6.(q) 



(2n 



{x ^ y)] , 
(A.45) 



where matrix G is defined as 
1 



G 



{,]' + J — J cos qyaz — 5J smqya2) (A. 46) 



and the determinant A equals 

A(qj,) = J'{J' + 27) + (J^ + {5Jf) sin^ qy . (A.47) 

When J = (5J = we find 

i^[b(q)] ^-^6a(q)&a(-q) , 

which restores the limit of an ordinary 2D XY model 
(r05ll . 

Returning to the partition function (|A.33I) and using 
the expression ljA.45p for i^[b] wc just found, we are now 
in position to integrate out the gauge field A(p) and ob- 
tain the analogue of the Coulomb gas representation for 
our model. Note that HA.30() implies 

6,(q) = (l-e-'^«)A(q) (A.48) 
6,(q) = -(l-e-^0A(q) . (A.49) 

The partition function now becomes 

/oo 
-oo p 

I H''("q)^(q)-A(-q)MA(q; 



i(p) 

exp 



(A.50) 



where 2x2 matrix M is given by 
1 



M 



2A(%) 



(1 - (73 cosgy)[j' + J(l - 0-3 cos Qy)] 



+ ai5 J sin"^ qy)^ +{x^y) . (A. 51) 
After integration over A(p) we obtain 



Z = ^ exp 
Hp) 



T 



2F(-q)M(q)[(q) , (A.52) 



where matrix M is the inverse of M . The elements of 
matrix M satisfy the following simple identities: 



A/ii(q) =M22(q-g) 
Mi2(q) =M2i(q) . 



(A.53) 
(A.54) 



Consequently, the integrand in the exponent of partition 
function IIA.521) can be written as 



(A.55) 



2 ;(-q)Mn(q)/(q) + ?(g - q)Afi2(q);(q) 



The explicit form of M — is rather cumbersome. 

Fortunately, we will only need the leading and subleading 
order terms in the long wavelength [q — > 0) expansion: 



Mn(q) 



+ 



4( J' + 2 J) 
{J' 



AJ+18SJ){qt+ql) + 12SJqlql 



iJ'{ql 
+ 0{q^) 
Mi2(q) - '5J + 0{q^) . 



(A.56) 
(A.57) 



The terms of order 0{q^) correspond to M{p — p') de- 
creasing at least as fast as \p — p'\~^. Returning to the 
real space representation, we obtain 

5^exp{-^ J2{1{P)MMP - pW)+ 



i(p) 



where 



p,p' 



^^■Pl{p)MMp-p')lip'))} , (A.58) 



mmp-p') = r ^ r ^e^-^-'^-^'hiM^) . 

(A.59) 

The two terms in the exponent of (|A.58p are easy to in- 
terpret. Recall that integers l{p) coupled to A(p) effec- 
tively reside at the centers of the plaquettes of the blue 
lattice corresponding to either black or red sites in Fig. 
121 The terms containing il/n clearly describe the av- 
erage interaction between two plaquettes irrespective of 
their "color" , while the terms with M12 reflect the differ- 
ence between the red and black sites. For example, the 
strength of interaction between two black sites separated 
by two lattice spacings with p — and p — 2x will be 
different from interaction between two red sites aX p — 
p' = ix due to the factor exp(ig • p) that multiplies Mi2. 

At large distances the Fourier transform can be evalu- 
ated by comparison to the standard lattice Green's func- 
tion in two dimensions. The difference 

il/ii(q) - 4(J' + 27)- ^- 

4 — 2 cos — 2 cos qy 

is finite at q = 0, and therefore the Fourier transform of 
this difference vanishes at large distances. Thus 



Mii(p) = 4(j'-F2J)x 

'"^ dqx_ r dqy 
2tt 



27r 4 — 2 cos qx — 2 cos qy 



(A.60) 
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Using the well known asymptotic behavior of the last 
integraS2>, we find 



Mn{p - p') ^ Mn{0)- 
_ A{J' + 2J) 
2tt 



[\n\p-p'\ + {) + ...] , (A.61) 



where Ci can be related to the Euler-Maschcroni con- 
stant 7 « 0.5772 as Ci = 7 + ln(2\/2). Note that 
Mii(p) formally logarithmically diverges because Mii(q) 
is proportional to at small momenta. The difference, 
Afii(p) — Mii(O), however, is finite. The overall infinite 
additive constant has a simple physical interpretation, 
just like for an ordinary two-dimensional XY model, as 
will become clear in a moment. 

The real space expression for M12 can be easily calcu- 
lated directly: 



Mi2{p- p') 



-5 J 60 



(A.62) 



where . . . denotes terms that decrease at least as fast 
as |p — p'l^^. Combining l|A.59l) . (|A.61|) . and ljA.62|) . we 
obtain the following expression, after separating off the 
terms with p = p': 



Z = ^exp[-^(Mn(0)^;2(p)- 



i{p) 



J2 M,,{p-p')l{p)l{p') - 5jY,e^^"l\p) 
p^p' P 

(A.63) 



By applying the long distance expansion of Mii(p) in 
the sum containing terms with p ^ p' we find that the 
partition function Z becomes 

^ exp [- ^ (Mn (0) ^(P)) ' " E ^'^"^^ 
i(p) P P 



4( J' + 2 J) 
2^ 



Y,l{p)l{p'){\n\p-p'\ + C, 



p^p' 



(A.64) 



We now return to the discussion of the formally di- 
vergent constant A/ii(0). This divergence is a reflection 
of the logarithmic dependence of a single vortex energy 
on the system size. If the number of the sites iV were 
finite, we would have obtained a constant of order IniV 
for Mii(O) instead of an outright divergence. Although 
finite, this constant becomes large in the thermodynamic 
limit iV — > c», with the effect of suppressing all configu- 
rations of the integer- valued field except those that 
satisfy 



(A.65) 



This is nothing but the charge neutrality condition in the 
partition function of a 2D Coulomb plasma. Restricting 



ourselves only to such configurations, we obtain: 
-7r2 /4(J' + 2J) 



Z = 



Hp) 



27r 



j2iip)Kp')Hp-p'\ 



Pt^P' 



+ + W^e'SPZ2(p)jJ (A.66) 
p 

A further simplification is achieved by noticing that 

i{p)Kp') ^ (T.^{p)) ~Y.^\p) ■ 

pi^p' \ p / p 

Since the first term on the right hand side vanishes by 
virtue of ljA.65|l . the partition function equals 

^ = E -P [ Y (^^^^^ E ^^P)W In \P P'\ 



Kp) 



Pt^P' 



4Ci(J' + 2J) 
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J2l\p) + 6jJ2e^^-^l\p) 
p p 

(A.67) 



This is the Coulomb gas representation of our model de- 
scribing charges l{p) residing on black and red plaque- 
ttes and interacting with long-ranged forces. Condition 
(| A. 6 5(1 therefore is simply an expression of the overall 
neutrality of the system - only the configurations with 
the same number of vortices and antivortices contribute 
to the partition function. 

The Hamiltonian of the system can be finally recast as 

nt = -n{j' + Ji + J2) ^ l{p)Kp') In |p - p'\ 

Pt^P' 

+ e:Y.I\p)+eIY.I\p) , (A.68) 

pen peB 

where the core energies of the vortices on red {TV) and 
black {B) plaquettes are expressed through the original 
parameters of the model as 



El = ^Ci{J' + Ji + J2) ~—{Ji- J2) (A.69) 



El = nCi{J' + Ji + J2) + —{Ji- J2) 



(A.70) 



The Hamiltonian HA.68() is of the form equivalent to Eq. 
(|23|l derived in the main text from the continuum formu- 
lation. Note that in the "low temperature" limit J ^ 1, 
coefficients J and J' = Pv{J) coincide (see Eq. (IA.21|) '). 
and the agreement with the continuum formulation is 
complete: the effective strength of the long range inter- 
action between vortices is 

J = J + J1 + J2 . 
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6. Quantum phase fluctuations 

The derivation in 2+1 dimensions follows closely the 
steps of the two-dimensional case considered in the pre- 
vious subsection. Denoting the imaginary time by r, the 
partition function of the model is 



where p is defined precisely like in the 2D case and r = 
{p,t). The Lagrangian of our quantum model is defined 

as 

-L{p,t) = +i7'/'r + j(cOs(Va:0r)+COs(Va0r)) 

i>r+y — (f'r+x))- 

(A.71) 



+ -^12(p)(cOs((/)r+x+ji - 0r) + ' 



The sign of the Berry phase is chosen to be positive for 
later convenience; obviously the partition function is not 
affected by the change. It is convenient to replace the 
integrals over continuous variable r by sums over discrete 
T„ (abbreviated often as r below) separated by intervals 
of "length" e. For brevity, we will use (pr+f to denote 

The terms containing time derivatives can be trans- 
formed as follows: 



exp 



(A.72) 



After completing the square we have 



exp^ 



2e 



Note that this expression can be formally replaced by a 
sum 



m{r) r 



2 ■ 



since, clearly, only the term m = survives in the limit 
of small e. The latter form is convenient because now the 
Poisson identity 



a 



-n + inc 



n— — oo 



E ^^P 



m— — oc 



27rr7i) 



2a 



(A.73) 



can be applied. The result is 

E ^^p[E 



exp 



^-E( ''''' 



2K, 



■Ur{vf 



(A.74) 



where Ur(r) is an integer- valued field. Using the identity 
the partition function Z can now be rewritten as 



EE / W_d4>r eyi^\i^{u^(Y)V 

u w± r r 

+ W+ (r) {(t>v+x+y - cf)v)+W^ (r) {(j)r+y - (j>v+x ) 

exp [-J] 



uI{y) , w\[v)+wl{r) e{ur~ff 



2J' 



2J{2(P) 



2A'o 



(A.75) 



The reader should bear in mind that throughout the ap- 
pendix the Greek indices exclusively denote the space- 
like components x, ?/ of a three- vector, while latin indices 
denote both space-like and time-like components, as the 
case may be. The coefficients J' and J'i2{p) a-re defined 
as 



J' = l3vieJ) 
Ji2{p) = Pv{Ji2{ep)) . 



(A.76) 
(A.77) 



We now proceed to transform the above expression by 
shifting the differences of the phases 4>r onto the dif- 
ference of the fields u(r) and w±(r) by using discrete 
integration by parts (jA.26p . just as it was done in the 
two-dimensional case: 

^ = E E / n • "(r) + 



exp 



w+ (r — i: — ?/)) + (w_ (r — x) 
ui{v) , eK-7)2 



-(r-y)) 



[-El 



2J' 



2Kn 



(A.78) 



where bold letters stand for three-dimensional vectors 
and V • u denotes three-dimensional divergence. 

In the absence of the next-nearest coupling terms rep- 
resented by w±, the 5-function constraint in (|A.78p is 
resolved by u = V x a, where the lattice curl is defined 
as 



eyfeVjAfe(r - Bk) 



In our case of a d-wave superconductor and finite 'w± , we 
rewrite the constraint as 



V • (u(r) + w[w+, w_] = 



(A.79) 
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where w = (wx,Wy,0) and Wr^(^y^ are defined in the pre- 
vious subsection section (IA.28|I . The solution is clearly 

u = V X A — w[w+, w_] . 

One can easily check that the choice of A is not unique: 
for arbitrary scalar function ^(r) 

[Vx (A + V^(r))]^- [Vx A]^ = 

= e^jfeVj VfeC(r - Sk) = e.jTc Vj- VfeC(r) = (A.80) 

This gauge invariance implies that a gauge-fixing term 
must be introduced when replacing the sums over integer 
field u by summation over A in order to avoid multiple- 
counting. The next step is most easily derived in the 
temporal gauge: 



<5(A3)^n^. 



A3(r),0 



Afterwards, the results will be generalized to an arbi- 
trary gauge-fixing condition. We proceed by rewriting 
the partition function as 



/(Vx A-F[w±])2 



A w± 



2J' 

,,2 



fA3(r),0 
(A.81) 



and apply the Poisson formula in order to obtain a the- 
ory depending on continuous rather than integer valued 
gauge field A: 

J2 <5A3,o/(Ai(r), A2(r), A3(r)) - / dAi(r)dA2(r)x 

A(r) 

J2 e2"S.;=(r)A„(r)j(A^(j.),A2(r),0) 
U(r),i2(r) 

/OO 
dAi(r)dA2(r)dA3(r)5(A3)x 

I- 

5^%,e2-£^'^W^^«/(Ai(r),A2(r),A3(r)) . (A.82) 

l(r) 

In performing the last step, we formally introduced /3A3 
and an additional sum over hir). The delta function 
(5(A3(r)) ensures that the exponent is not affected. All 
terms in the sum over I3 are therefore equal, and in or- 
der to avoid multiple-counting we need to impose a con- 
straint, chosen as V • Z = 0, by assigning 

^3 = -(V3)"\Vi/i+V2?2) . 

Note that the result of applying operator (V3)~^ to an in- 
teger field is another integer. The integer-valued field l(r) 
with zero divergence describes non-backtracking closed 



loops on the 2-1-1 space-time lattice. The field A(r) is 
now continuous and the temporal gauge condition 



e.[A]=n'5(A3(r)) 



can be replaced'^^ by an arbitrary gauge-fixing condition 
0[A]; examples are V • A = (Landau gauge) or • 
A^ = (radiation gauge): 

^JA3,o/(Ai(r),A2(r),A3(r)) = 

A(r) 

00 

nrfA(r)e[A(r)]^e2-5:.Kr)-A(r)^(^(^))^_^ _ 

r l(r) 

(A.83) 

This final identity allows us to rewrite our partition func- 
tion as 



re 



i:^F[(VxA)j 



J|dA(r)e[A(r,^ 

J —00 J. 

^_^gi:42"'(r)-A(r)-(e/2ifo)((VxA)o-7)'] 



l(r) 



where F[{W x A)±] has already been calculated in ljA.45|) 
and can be used as is, provided that proper definitions 
(|A.77I) of J' and Ji2(p) are replaced. 

The remaining steps leading to the "Coulomb" gas rep- 
resentation of 3D vortex loops are conceptually similar to 
the 2D case from the previous subsection. The algebra, 
however, is considerably more involved. We therefore 
will go slowly and first wade through the derivation for 
the simple case Ji = J2 = 0. This is just the ordinary 
(2-|-l)D XY model, appropriate for our s-wave pedagog- 
ical exercise from the main text and the beginning of 
this Appendix (|A.2p . Only the configurations w±{r) = 
contribute to the functional F[(V x A)±] l|A.45l) and we 
recover the usuali^i anisotropic 3D XY model in a uni- 
form magnetic field H = /f: 



Y[dAir)e[A{r)]J2S- 



v-i- 



l(r) 



[27r^i(r)-A(r)- <^|^-5^ ((VxA)o-/) 



(A.85) 



To obtain the lattice loop gas representation, we need 
to integrate out the gauge field A. The most transparent 
connection with the results for the 2D model is obtained 
by using the radiation gauge VqAq, = 0: 

Y,{WxA)l^Y.^W.Ay{r-y)-WyA.4r~x)f (A.86) 

r r 

Expanding the square and shifting the difference opera- 
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tors via IIA.2611 we have 



^(V xA)l = J2 -^yi^ - y)V,V,Aj,(r - y) 

r r 

- A^(r - x)'SIy^ykx{v -x)- TsJyKy{v - y)'^xKx{T - x) 
= ^-A^(r)V,V,A^(r)-(V^A^(r))2 . (A.87) 



We introduce the following notation for the lattice ana- 
logues of wavevectors qj appearing from discrete left- 
sided or right-sided derivatives after Fourier transforma- 
tion: 



Q.(q) = 



e»9j - 1 



1 - e-"i^ 



Qf(q) = Q,(q-g) 



— E — 1 + e"*«J 



(A.88) 

(A.89) 

(A.90) 

(A.91) 
(A.92) 



The arguments of Qj and is assumed to be q unless 
specified otherwise. 

We define the Fourier transformation as 

where the sum over frequencies go runs through go ~ 
0, . . ., ^ and the integrals over qx^qy extend from — tt 
to TT. Using the definition and properties that follow from 
it 

/(q) = ^E / ^^"'""/(q) (A.93) 

we obtain in the radiation gauge: 
^(V X A)2 = - ^A^(r)V.V.A^(r) 

r r 



(A.94) 



= ^3J)2^M(-q)Q.Q.A4q) , (A.95) 

Similarly, 

^(V X A)i = ^(-Ao(r)V,V^Ao(r) 

r r 

- A^(r)VoVoA^(r) - 2VoAo(r - f)V„A„(r - e„) ) , 



(A.96) 



where the last term vanishes due to our choice of the 
gauge. In the momentum space we have: 

$:(VxA)i^^ j ^|3i_(A„(_q)Ao(q)Q,Q, 



A^(-q)A^(q)QoQo • (A.97) 



The above definitions are generally valid but now we 
focus again on the simple case of Ji = J2 = 0. The 
partition function 2'o, given by l)A.85(l . can be written as 



E h-i^ 



l(r) 



n^Ac 



r X 



(5[VaAJe 



i:^[2.»i(r).A(r)-<I|i^-5^(VxA)^J 



(A.98) 



In arriving at the expression above we performed a shift 
of variable A — > A -|- Ay, where Ay is a time independent 
vector potential corresponding to a constant and uniform 
magnetic field /f. After Fourier transformation the ex- 
pression in the exponent can be written as 



E 



dqj 



/3e(27r)2 



27rj;(q) • A(-q) - 
1 „ ^ 1 



Ao(-q)Ao(q)Q^Q^ 



2J' 



A^(-q)A^(q)(— QoQo + ^Q.Q.)] ■ (A.99) 



Note that temporal and space-like components are inde- 
pendent and can be integrated out separately. Integra- 
tion over Aq is trivial and yields 



exp 



^2vrV'E 



«0 



dqxdqy Zo(-q)'o(q) 



exp 



[-2.V'^ 



dqxdqy /o(-q)'o(q) 



/3e(27r)2 4 — 2 cos — 2 cos qy 



(A.lOO) 



The remaining integral 

n dA(q)<5[Q„A„] exp (2^*;,(-q)A,(q) 



[jJjQoQo + ^Q.Q.)A^(-q)A^(q; 



(A.lOl) 



can be computed by switching to 2D transverse and lon- 
gitudinal components of A, which we define on the lattice 
as 



Ai(q) 



AT(q) = i 



^Q,Ax{<l)+QyAy{q) 

Q± 

. -Qy Ax{q) + QxAy{q} 



Qa 



(A.102) 
(A.103) 



where 



Q± = yQc.Qc.^ jQxQcc + QyQ 
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The Aa; (r) and Ay (r) can be expressed through the trans- 
verse and longitudinal components of the gauge filed A 
as 



A.(q) 



Aj,(q) = -i 



. QxAL(q) -Q^AT(q) 
QyAi(q) + g^AT(q) 



(A.104) 



(A.105) 



Now observe that the Jacobian of the transformation 
(Aa;, Ay) (Ai, At) is unity, 2D divergence of A is pro- 
portional to Ai as expected: 



and 



;„(-q)A„(q) = ZT(-q)AT(q) + ;L(q)Ai(q) . 
The integral IjA.lOip can be written as 

/_^n^A.(q)exp[;^ / ^(2..Z,(-q)A,(q) 



^QoQo + l^Qo^Qc ) AT(-q)AT(q) 



(A.106) 



which finally gives 
Zo oc ^ exp 

l(r) 



/ aqxdqy f 
^ ^.1 e/3(27r)2V 



dqxdqy ( ^jiki-aijlolq) 



^T(-q)^T(q) 



QaQa 

. (A.107) 



2J' QoQo ~l" 2K' QaQa ' 

This is the desired vortex loop gas representation of our 
model. Rather than integrating out the gauge field A in 
(|A.85p . one can partially perform the sum over over the 
integers l(r) and arrive at yet another (dual) representa- 
tion of partition function Z . The constraint V-l(r) = in 
the partition function Zq (jA.85|l is rewritten using aux- 
iliary variables a(r) as 

+ 27ri l(r) • A(r) + i ^ a(r)V • l(r) 

r 



l(r) 



n 



+ i 



^l,(2^A,(r)- V,a(r)) 



27r 



da{r) exp 



^/3'cos(V,a(r)-27rA,) 

(A.108) 



The last equation describes a lattice superconductor. 
Note that the partition function (|A.85p does not contain 
quadratic terms l^(r). Instead, one introduces a vortex 
core energy term 

I Vl'(r) 

2/3y(/3') ^ ^ > 

by hand, and then, in the final lattice superconductor 
representation, a limit /3 ^ oo is taken. Such a system is 
called a frozen superconductor. Alternatively, terms pro- 
portional to P can be kept finite. Such "unfrozen lattice 
superconductor" is equivalent to an ensemble of vortex 
loops, namely an XY-model augmented with an addi- 
tional core energy that makes formation of vortices more 
difficult. Thus, applying HA.108|I to HA.85|) we obtain the 
partition function Zq describing lattice superconductor 
in a field / coupled to fiuctuating gauge field A: 



n 



da{r) / dA(r)e[A(r)]x 



exp|^^/3'cos(ViQ;(r) - 27rAi) 
5:f^(VxA)i + ;^((VxA)o-7)^ 



(A.109) 



The last step of our derivation is the standard^^ 
Ginzburg-Landau expansion of the action and for com- 
pleteness we reproduce here the derivation following 
Kleinert^^. 

First, we introduce a complex field — exp(iQ;(r)) 
and define covariant derivative operators Di, Di accord- 
ing to 



^^^^(r) = $(r -I- i)e-2"^- - $(r) 
D^^{r) ^ $(r) - $(r - x)e2'^*^-('-*) 



(A.llO) 
(A.lll) 



The following identity, which expresses the cosine in 
(jA.lOQp through Ur, can be proved easily: 

^cos(V,a(r) - 2^A^a;) = 5I^r(l + ^DMUr ■ 



Thus, for a given fixed configuration of dual gauge field 
A in (|A.109|I . the sum over all configurations of angular 
variables a(r) is 



cos(ViQ(r) — 27rAi) 



Zxy[A] = J PaWe^-''' 
and can be transformed into 

Zxy[A.]^ f VaexpiSP'^U^VUr] . (A.112) 

where V = (l + ^DiDi). Operator V is Hermitian, and 
therefore allows decomposition V = K? . Let us show 
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that ZxY is proportional to 

(A.113) 



where notation 

is used. To establish the equivalence, we group the terms 
in the exponent as 

- ^($; -6/3'5^C/;Xr'r)(<i>r " 6/3' /Cr'rC/r' ) 

J.! J./ 

+ 3/3'^^C/;„/Cr"r/Crr'f/r' ■ (A.114) 
r r'r" 

After a shift of variables and integrating out the auxiliary 
fields $r the result coincides with IIA.112|l up to an unim- 
portant proportionality factor. To obtain the equivalent 
description in terms of field we now integrate over 
the angular variables a-r ■ To simplify notation, we define 
Xi = Ai<f>/2 and X2 = K.^^* /2, or more explicitly. 



The partition function IjA.llSp is the desired dual rep- 
resentation of our initial anisotropic XY model with the 
Berry phase, in the simple case Ji = J2 ~ 0. 

Armed with the experience from the above derivation 
we now return to the partition function Z of the full- 
fledged model ljA.84p containing Ji and J2. As we will 
demonstrate, the effect of the next nearest neighbor in- 
teractions will be rather modest: to the leading order 
only the term proportional to |$(r)p will be modified. 
The prefactors of this term will be modulated, having 
different values on the red and the black plaquettes. 

To arrive at the dual representation of Z, we seek a 
gauge that will ensure the decoupling of the temporal 
and spatial components of A(q), similarly to the radi- 
ation gauge in the simple example above. The bilinear 
terms in A^ appearing in the exponent can be classified 
as following: first, there is a contribution from the term 
(V X A)q which in an arbitrary (yet unknown) gauge has 
been already calculated in (|A.87|I 

-^(A^(-q)Ap(q)Q,Q. - Q,Q^A,,(-q)A,(q; 

(A.116) 

To facilitate the bookkeeping of various terms resulting 
from F[(V x A)_l] we use the following set of identities: 



QoA^(~q)CT3 
(A.117) 



Integration over the phases in ljA.112|l now amounts to 
calculation of disentangled integrals at separate r: 



where /g denotes the modified Bessel function. Thus, 
omitting non-essential overall prefactors the expression 
for ZxY [A] assumes the following form: 



Finally, after applying the Taylor expansion 



2 4 



and retaining only the leading terms, we obtain 
Zo= f [|d$(r)d$*(r)dA(r)e[A(r)]x 
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2K' 







Ao(q)-ie "^"Qo^sA^lq) 



(A.118) 

The expressions for fey(q) can be obtained by replacing 
X 'r-t y and the overall change of sign. Using these iden- 
tities the bilinear form 



) (Gniqy) ^ ( ) + (a; ^ y) 



can be written as a sum of two groups: the diagonal 
terms are 



}T, I Q_yQyGii{qy) QyQyGl2iqy)\ , 

' \QtQIG2l{qy) QtQIG22{qy)J 

+ A'^(-q)\ *5oQoGii(gy) ~_Q oQnG 12 (qy) \ , ^ 

' \-QoQoG2l{qy) QoQoG22{qy) J ^ 

+ {x^y) . (A.119) 

In addition, we obtain cross-terms that couple the spatial 
and temporal components of A^ 



(A.115) QoA'^{-q)(P{qy)Ay{q) + P{q^)AA<l))+c.c. , (A.120) 
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where c.c. denotes complex conjugation and matrices 
P{qa) are defined as 



QlG2l{qa) -QlG22{qa)j 



Note that the terms diagonal in Aq are exactly what we 
encountered in (jA.51|) when we considered the 2D exam- 
ple. The ofT-diagonal terms can be eliminated altogether 
by choosing a gauge defined by the following relation be- 
tween Aj;(q) and Aj,(q): 



P{qy)Ay{q) + P{q,)kAq) = 



(A.121) 



The matrix equation can be resolved by Aj,(q) = 
r(q)Aa;(q) where matrix F is defined as r(q) — 
—P~^{qy)P{qx)- The spatial part of the action in mo- 
mentum space becomes 



E 

10 



dqxdqy 
/3e(27r)2 



(/I(-q) + [^(-q)r(q))A.(q) 

-A,(-q)l(q)A,(q)l , (A.122) 



where 2x2 matrix X is defined as 



1 





r^(-q 



Gii{qx) -Gi2{qx) 
-G2l{qx) G22{qx) , 



Gii{qx) -Guiqx) 

-G2l{qx) G22{qx) 



r(q) . (A.123) 



Note that we omitted the second term in (lA. 116(1 as it 
is only of the order Q\ and can be safely neglected for 
extracting the long distance behavior. After integrating 
out A^;, we obtain 

e^P[ 4 2^7 /3e(27r)2'' 
go 

(/I(-q) + f^(-q)r(q))x(4(q) + r^(-q)[,(q))] , 

(A.124) 

where X — X^^. Expanding this expression in the region 
of small momenta as in (|A.55p and retaining only the 
leading order terms, as in the two-dimensional example, 
we find 



_ 2V^ f dqxdqy 1 

" H2n)^ijll+2h^r^'' 

(qylxi-q) - qxlyhnij (qylxiq) ~ qxlyiq) 



(A.125) 



Observe that there are no cross-terms that couple modes 
at wavevectors q and q — g to the order of and q^^. 
The spatial part of the action after integrating out the 
gauge fields is equivalent to the result l|A.107|l obtained 
in the framework of simple model Zq where J' is replaced 
by the effective coupling constant J = J' + Ji + J2 • 
Combining ((A. 125(1 and ijSI 



the final form of the 
action Z in terms of closed vortex loops l(r) can be writ- 
ten as 

Z = ^%.,exp[^(2«l(r) • A/(r) - E',{p)llir) 

l(r) 

2Y^ f dqxdqy ,~lo{-q)lo{cL) 



90 



/3e(27r)2 



9i 



ZT(q)/T(-q)) , (A.126) 



where E'^{p) = ±7r^(Ji — J2)/4 depending on whether p 
corresponds to a black or a red plaquette. 

We had intentionally used momentum representation 
for the last two terms in the exponent. It is important to 
recognize that these terms are precisely what one would 
have obtained for the usual 3D XY model with no next 
nearest neighbors interaction and the effective nearest 
neighbors coupling constant equal to J — J' + Ji + J2- 
Thus, we may introduce a dual gauge field A(r) and 
present the partition function as 



Z 



E 

l(r) 



nrfA(r)e[A(r)]%.,x 

00 J. 

exp[^(2^* l(r).(A(r)+A/(r)) 



By shifting A ^ 



(V X A)l (V X A)^ 



2J 2K' 
A — A/ back we obtain 



(A.127) 



Z = 



Y[ dA(r)e[A(r)] exp (2Tri l(r) • A(r) 

r r 

(VxA)i [(VxA)o-7]2 



2J 



2K' 



(A.128) 



Note that in the absence of the next nearest neighbors 
interactions, a similar expression l(A.85() contained no 
quadratic terms P(r) and the core energy — ■g^jj^l^(r) 
was introduces by hand. In the present case, the differ- 
ence Ai?^ of the core energies on the black and red sites 
is finite due to the anisotropic next nearest neighbors in- 
teractions. However, the average magnitude is still zero 
within our model, and here we also need to introduce a 
constant average core energy term ^^^^,^ P(r). Thereby 

we replace E'^{p) in lA. 128(1 by 
1 1 



+ K(.P) 



41 



where the function (3'{p) is implicitly defined by this 
equation. 

The remaining steps repeat the derivation leading from 
(|A.108|) to HA.115|) with the replacement — > /3'(p) and 
result in the Ginzburg-Landau expansion of our dual the- 
ory: 

Z = y"]Jd$(r)d$*(r)dA(r)e[A(r)]x 



exp 



This is our final result - the partition function (jA.129|) 
represents the Ginzburg-Landau functional of a dual 
type-II superconductor appropriate for our model and 
subjected to a constant dual magnetic field /. 



1 f 1 



l)l'i>Wl' + ^l<I'(r)|^ 



^(VxA)i + ^((VxA)o -/)'}] . (A.129) 
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